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ABSTRACT 

We present the highest resolution, wide-field radio survey of a nearby face-on star¬ 
forming galaxy to date. The multi-phase centre technique is used to survey the entire 
disk of M51 (77 square arc minutes) at a maximum resolution of 5 milli-arcseconds 
on a single 8 hr pointing with the European VLBI Network at 18 cm. In total, 7 
billion pixels were imaged using 192 phase centres that resulted in the detection of six 
sources: the Seyfert nucleus, the supernova SN 2011dh, and four background AGNs. 
Using the wealth of archival data available in the radio (MERLIN and the VLA), 
optical {Hubble Space Telescope) and X-rays {Chandra) the properties of the individual 
sources were investigated in detail. The combined multi-wavelength observations reveal 
a very complex and puzzling core region that includes a low-luminosity parsec scale 
core-jet structure typical of AGNs, with evidence for a lateral shift corresponding to 
0.27c. Furthermore, there is evidence for a fossil radio hotspot located 1.44 kpc from 
the Seyfert nucleus that may have resulted from a previous ejection cycle. Our study 
provides measures of the supernova and star-formation rates that are comparable to 
independent studies at other wavelengths, and places further limits on the radio and 
X-ray luminosity evolution of the supernovae SN 19941, SN 2005cs and SN 201 Idh. 
The radio images of background AGN reveal complex morphologies that are indicative 
of powerful radio galaxies, and confirmed via the X-ray and optical properties. 

Key words: galaxies: Seyfert - galaxies: individual (M51) - techniques: radio astron¬ 
omy - interferometric 


1 INTRODUCTION 

The formation of stars is a fundamental astrophysical pro¬ 
cess that takes place in almost all galaxies along the Hub¬ 
ble sequence and as such plays a major role in galaxy 
evolution. Star formation can occur in isolated areas or 
over the entire galaxy, resulting in 6-7 orders of magni¬ 
tude variation in galaxy star formation rates llKennicuttI 
Il998ll . Galaxies that are continuously forming stars over 
an entire galaxy (hereafter star-forming galaxies) typically 
host rich populations of young stars and star clusters, HII 
regions, giant molecular clouds, and supernova remnants 
(SNRs). Furthermore, it is not uncommon for star-forming 
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gala xies to host low luminosity Seyfert AGN at their cores 
(e.g. iFord et al. 1^8^ Pellegrini et al. 2000j^ Lenc fc Tingayl 

200^ Perez;Torreset_alJ_ 201^ Rampadarath et al.l 201ol : 

Paggi et al.l l2013l : IPanessa fc Girolettil l2013ll . suggesting 
possible connections between star formati on and low lu mi- 
nosity AGN activity in these galaxies le.g. lNetzerll2009l l. 

Nearby star-forming galaxies are ideal laboratories 
for understanding the nature of supernovae (SNe), SNRs, 
and the nuclear environment on parsec and sub-parsec 
scales through the use of wide-field, high resolution ra¬ 
dio observations. With multi-epoch radio monitoring, it 
is possible to construct the star formation and SNe his¬ 
tory of these galaxies. Th is was attempted for the star- 
burst galaxy NGC 253 by iRampadarath et al.l (120141 '). us¬ 


ing data spanning 21 years. Similar studies have been 
conducted of nearby edge-on, star-forming galaxies: M82 
ll Pedlar et al.l 1 19991 : iMcDonald et al.1 12002 : iBeswick et al.l 
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200^ Fenech_et_^ l2010ll: Arp 220 (|Rovilos_e^^ 


Arp 299 du vestadi 


Lonsdale et al. 200d Parraet^ 2007 |jBateja^et^ 


I 2 OO 9 I: iR.oniero-Canizales et al. 20^; 

2 OI 2 I): NGC 4945 jlLenc fc Tingavl 1200911: 


■200d: 

l2012ll: 


iBondi et al 

NGC 253 


(Ulvestad k. Antonuccil 19971: Tingavl 2004 


iLenc fc TingavIboOfl l: and M31 llMorgan et alJ 201L ), and 


nearby face-on, spiral galaxies: M 83 ([Tamer_^Ho 199j; 
ICowair et al. 19941 : Maddox et al. 2006r Long_^£_^ 201 


M51 ll C rane fc van der HulstI Il992l: [Tur ner fc He 1 19941 : 
iMaddox et alJl2007l): NGC 300 I Pannuti et al. 200oli: NGC 
7793 ( Pannuti et alJ I 2002 I I: and NGC 6946 ( Lacey et alJ 
[1993). Furthermore, there have been numerous surveys of 
face-on, spiral galaxies with VLBI, focussiuK on the nuclear 
regio ns (e.R. Bietenlmlz_et^l.l 2004 ^ Giroletti_^Panesd 


2009|;__ |Borhempi et al.l l2012l : IPanessa &; Girolettil l2013l : 


Dpi et al. 20131 ). 


Here, the results of the highest resolution, wide-held 
radio survey of a nearby, grand design spiral galaxy are pre¬ 
sented. The target for this survey is the Whirlpool galaxy, 
(aka M51a, NGC 5194), that is currently undergoing an in¬ 
teraction with a smaller companion, M51b (aka NGC 5 195). 
At a distance of 8.4 ± 0.7 Mpe dVinko et al.l I2OI2I) and 
with a nearly face-on orientation, M51 is ideal for popu¬ 
lation and morphology studies at all wavelengths. Like most 
spiral galaxies M51a is actively forming stars predominantly 
(but not exclusively) within it s spiral arms dCalzetti et al.l 
I2OO5I : iKaleida fc Scow^ |201(]|) . with a slightly enhanced 
formation of young^ mas sive star clusters towards M51b 
dKaleida fc ScowenTuOl^ . Furthermore, in the last 70 
years M51a and M51b have ho sted four optically ob - 
served SNe: the Type la SN 1945A dKowal fc SargentHl97ll) : 


the Type Ib/c SN 19941 ( Puckett et al. M 19941 ): the Type 


Ib/c SN 20 05cs d Muendl ein et al.ll2005l) : and the Type Ilb 
SN 2011dh dCriga et al.ll201ll ). In addition, high resolution 
radio observations of the nuclear region with the VLA re¬ 
veals the presence of a bidirectional jet, which is associated 
with a low-luminosity Seyfert 2 AGN, located within the 
nucleus of M51 a, that interacts with the surrounding inter¬ 
stellar medium dCrane fc van der Hulstll 19921 : iBradlev et al.l 
l2004lMaddox et al.ll20(3^ 

The first high-resolutio n radi o survey of M51a was 
performed by iMaddox et al.l d2007l) at 20 and 6 cm, with 
the VL A. With a reso l ution of 1.5” x 1.2" at both wave¬ 
lengths, IMaddox et al.l d2007l ) detected 107 compact radio 
sources, distributed throughout the disk of M51a but pre- 
domi nantly concent r ated within the inner 2 kpc nuclear re¬ 
gion. [M^d^^et^^ (I 2 OO 7 I ) compared the results of the radio 
observations with multi-wavelength observations and found 
of the 107 radio sources: 44 radio sources with large HIT 
counterparts; 24 with stellar cluster counterparts; 13 with 
X-ray counterparts, most probably X-ray binaries or micro- 
blazars; and six radio sources associated with Ha emission, 
possibly young SNRs with resolved shells a nd the remain¬ 
ing as either background AGNs or unknown. IMaddox et al.l 
(l2007t) estimated the age for the SNRs to be ~ 2000-3300 


yr. 

Through application of the technique of wide-held 
VLBI, with the European VLBI Network (EVN), we sur¬ 
veyed the full disk of M51 (11'.2 x 6'.9) for faint, compact 
radio sources at 18 cm. Section[2]describes the observations, 
wide-held correlation, calibration and data reduction meth¬ 
ods that were employed within this study. In Section [31 the 
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Figure 1. Optical image of M51 from the Space Telescope Sci¬ 
ence Institute digitised plate archives (image credit: DS9 image 
server), overlaid with the sensitivity contours of the EVN. The 
contours are 6.7 pjy, 7.5 pjy and 13 pjy, and are the theoretical 
sensitivity taking into account the total observing time, band¬ 
width, the SEFD, effective area, data rate and the primary beam 
of the individual antennas. The red crosses are the positions of 192 
simultaneous multi-phase centres obtained through DiFX used to 
survey M51 with the EVN. 



wide-held VLBI imaging, source detection methods, and pri¬ 
mary beam correction are disenssed. Radio images of the 
sources detected are also presented. The multi-wavelength 
and multi-resolution ancillary data (radio, x-rays and opti¬ 
cal) used in this study are presented, including data reduc¬ 
tion steps, in Section [T] Discussion of the sources detected 
are given in Section[5] including comparisons with the ancil¬ 
lary data and previous results obtained from the literature. 
Finally, the study is concluded in Section (31 For power law 
spectra we adopt the convention that S . 

2 OBSERVATIONS, CORRELATION AND 
DATA REDUCTION 

2.1 Observations 

M51 was observed on 2011-11-7 for 8 hrs (UTC: 01:30 - 
09:30) at a wavelength of 18 cm with a single target point¬ 
ing of the EVN. The experiment was performed using the 
Effelsberg (Ef), Jodrell Bank 76 m Lovell dish (Jbl), a single 
Westerbork dish (Wb), Medicina (Me), Onsala 25 m (On), 
Svetloe (Sv), Zelenchukskaya (Zc), Badary (Bd), Urumqi 
(Ur), and Shanghai (Sh) antennas. During the observation, 
5 minute scans of M51 (centred on: RA = 13^29™ 52.698''; 
Dec = -|-47°ll' 42".930 [J2000.0]) were scheduled, alternat¬ 
ing with 1 minute scans of a phase reference calibration 
source, J1332+4722 (RA = 13'’ 32“ 45.246424"; Dec = +47° 
22' 22".667700 [J2000.0]) located 31' from M51. 

The VLBI data were recorded in 8x 16 MHz sub-bands, 
dual polarisation with 2-bit data sampling at the Nyquist 
rate, resulting in a total data rate of 1024 Mbps, producing 
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Figure 2. The (u, v) coverage of the 10-station EVN observa¬ 
tion at 18 cm (One point per sub-band per target scan, for each 
baseline only). 

a Icr sensitivity (at the pointing centre) of 6pJy. Figure [T] 
shows an optical image of M51 overlaid with the sensitivity 
contours of the EVN observation, taking into account the 
different primary beam of each antenna. The angular size of 
M51 (11'.2 X 6 '.9) allows most of the galaxy to lie within 
the primary beam of the 100 m Effelsberg antenna (~ 7'.5). 
Eigure[2]shows the (u,v) coverage of the EVN observation. 
The maximum projected baseline is approximately 45 MA 
(Jbl-Sh), while the shortest projected baseline (Ef-Wb) is 
1.5 MA, hence structures with an angular size larger than 
150 mas (5.7 pc at the distance of M51) are resolved-out 
(i.e. structures with intensity variations on larger angular 
scales are not detected). 

2.2 Correlation 

To survey the entire disk of M51 for radio sources with the 
full capabilities at VLBI resolution requires a pplication of 
the t e chnique of multi-pha ae centre correlation llPeller et al.l 
I 2 OIII : [Morgan et al.ll201ll ). This was achieved through the 
use o f the DiFX software correlator llDeller et al.l l2007l . 
[ 2 ^, at the Max-Planck Institue for Radio Astronomy, 
Bonn, Germany, as at the time of this experiment the 
JIVE software correlator (SEXC) was not offering multi¬ 
phase centre correlatioiQ. Prior to correlation, the data were 
shipped via courier from each station to Bonn, where the 
correlation was performed with a high performance com¬ 
puter cluster comprising 60 computer nodes (8 cores each). 

The correlation was initially performed with sufficient 
frequency and time resolution (488 Hz and 1 ms) to restrict 
amplitude losses below 5% at a radius of 6' on the longest 
baseline due to bandwidth and time smearing). Then using 
the simultaneous multi-phase correlation method, the field 
was separated into 192 individual phase centres, allowing 


the area around each phase center to be imaged separately. 
The locations of the phase centres are shown in Figure [T] 
To minimise cumulative amplitude errors to less than 
5% due to non-coplanar baselines, bandwidth, and time- 
average smearing on the longest baseline per phase centre, 
each sub-band (of 16 MHz bandwidth) was averaged to 128 
channels (Si/ = 125 kHz) and time resolution of 1 second, 
resulting in a data size of 41 GB per phase centre. The cor¬ 
related data were transferred to the iVEC Petabyte Data 
Stor^l located at Curtin University, Perth, Western Aus¬ 
tralia, using the network data transfer protocol gridFTlQ- 


2.3 Data Reduction 

2.3.1 Data Preparation 

During correlation, the clock at the Badary antenna expe¬ 
rienced a glitch that caused it to skip ahead by Is. This re¬ 
sulted in a separation of the final correlated data per phase 
centre, at the time of the clock skip, into two individual FITS 
files. Concatenation of the two correlated files per phase cen¬ 
tre to form the individual visibility FITS via the aips task 
FITLD takes ~2 hours on a standard iCore 7 desktop com¬ 
puter. Extending this process to 192 phase centres would 
require more than 300 hrs of computing time, which is al¬ 
most tripled when considering read and write time of the 
FITS data to disk. 

In order to concatenate, calibrate and image the 192 
phase centres in a time that is suitable for this project, the 
Curtin ATNF VLBI Experiment (cave), which is a high per¬ 
formance computing cluster located at the Australia Tele¬ 
scope Compact Array (ATCA), Narrabri, NSW was used. 
CAVE comprises 14 computer nodes, where the nodes are 
quad-core 2.7 GHz AMD Sunfire x86 machines, with 16 GB 
memory and 400 GB disk capacity each. The main disk space 
(8 TB X 3) is located on a single centralised node (here¬ 
after, cave-store). The cluster is owned by Curtin Univer¬ 
sity and operated by the Australia Telescope National Fa¬ 
cility (ATNF). To use the multiple nodes of cave effectively 
for the concatenatio n, we made use of th e ParselTongue 
scripting language llKettenis et al.|[2006l) . ParselTongue 
allows the execution of multiple aips tasks simultaneously 
on a number of remote computers via the paralleltask 
module. 

The correlated datasets were copied from the iVEC 
Petabyte Data Store to the 8 TB data storage on cave- 
store. Before starting the concatenation process, the 
ParselTongue server is started on the individual node^. 
Executing the ParselTongue script on cave-store be¬ 
gins by copying the datasets to the individual nodes, such 
that the datasets for the N*** phase centre is sent to the N*** 
node. After the data are copied, fitld loads and concate¬ 
nates both datasets to aips and the result written to disk 
(via FITTP), copied to cave-store and eventually is copied 
to iVEG. Note, the data processing on the individual nodes 
is executed in parallel. 


® http://www.ivec.org/ 

^ The SFXC is currently capable of correlating multiple phase ^ http://www.globus.org/toolkit/docs/latest-stable/gridftp/ 

centres and this is offered as a standard EVN mode. ® See the ParselTongue online documentation 

^ http://www.atnf.csiro.au/people/Emil.Lenc/Calculators/wfcalcIpthgi://www.jive.nl/jivewiki/doku.php 
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Due to the limited storage space and network connec¬ 
tion between iVEC and cave-store, the datasets were sep¬ 
arated into ~1 TB groups (~ 25-30 phase centres) for pro¬ 
cessing. The time taken to process ~1 TB of data on two 
nodes on cave was 14 hrf0, resulting in a total time of ~87 
hrs to concatenate the 192 phase centres. 


2.3.2 Calibration 

The data reduction and calibration were performed using 
the data reduction package aips. Prior to calibration, flag¬ 
ging of data during times at which each of the antennas were 
known to be slewing and time ranges that contained known 
bad data was carried out via application of flag files and 
information provided in the observing logs. Corrections to 
the amplitudes in cross-correlation spectra from the auto¬ 
correlations were obtained via ACCOR. The amplitude cali¬ 
bration was obtained through application of antenna system 
temperatures (in Kelvin) measured during the observations, 
along with the gain (in Janskys per Kelvin) for each an¬ 
tenna. Following the amplitude calibration, a single value 
of the delay (i.e. phase against frequency) was determined 
for each antenna, polarization, and 16 MHz sub-band us¬ 
ing two minutes of data on the bright radio source 3C 345 
(aips task fring). Global fringe-fitting solutions were deter¬ 
mined for J1332+4722 (aips task fring) with a one minute 
solution interval, finding independent solutions for each of 
the 16 MHz sub-bands. The delay and phase solutions were 
examined and, following editing of bad solutions, applied 
to J1332+4722. The J1332+4722 data were then exported 
to DIFMAP (IShepherdI Il997fl . where the data were vector- 
averaged over 30 s, flagged of bad data, and imaged us¬ 
ing standard imaging techniques (deconvolution and self¬ 
calibration of both phase and amplitude). The resulting im¬ 
age of J1332+4722 is shown in Figure [3l The image shows 
no evidence for extended structure on these baselines at this 
frequency to a dynamic range of 2500, for regions in the 
image away from the inner side-lobes around the source. 
The final calibration solutions (phase and amplitude) from 
J1332-I- 4722 were exported vi a the difmap extension, COR- 
dumfQ dLenc fc Tingavl l2009l l to an AlPS-compatible solu¬ 
tion table. The solution table was then transferred to aips 
and applied to all sources in the dataset. The final calibrated 
J1332-I-4722 data were also used to derive a bandpass cali¬ 
bration via the aips task bpass. The edge channels of each 
band were flagged (5 channels from both the lower and up¬ 
per edge of each 128 channel band). The final calibration 
solutions were then applied for a first imaging pass for each 
phase centre. 


® The processing time per 1 TB group can be speeded up by 
using multiple nodes for the data processing i.e. three nodes would 
reduce the time to 9 hrs, four nodes to 7 hrs etc. However, cave, 
is a shared cluster and only two nodes were fully available for this 
project at any time. 

^ http://www.atnf.csiro.au/people/Emil.Lenc/tools/Tools/C 



Right Ascension (J2000) 

Figure 3. Contour map of the calibrator J1332+4722, with a 
dynamic range of ~ 2500 between peak flux density and the 
RMS for regions in the image away from the inner side-lobes 
around the source. The contours begin at ±1% of the peak 
(226.02 mjy/beam), and increase in increments of 2"% of the 
peak where n are integers (1...7). 



Figure 4. The distribution of the peak signal to noise ratio 
(S/N) for the 192 phase centres, where the x-axis is plotted in 
logscale. 

3 IMAGING, SOURCE DETECTION AND 

RESULTS 

3.1 Wide-field VLBI Imaging 

Prior to imaging, data affected by RFI were manually 
flagged via the aips tasks spflg and editr. Data received 
by the Me antenna were found to be corrupted by RFI. Ini¬ 
tial attempts at selective RFI removal were found to be ex¬ 
cessively time consuming and the entire antenna was flagged. 
Improving the final sensitivity by ~ 10%. Furthermore, data 
with u coordinate < 50 kA which are susceptible to RFI due 
to the low fringe rate were also flagged. The flag tables along 
with the final calibration and bandpass tables were copied 
and applied to all phase centres via split. 

Initial imaging was restricted to the four European an¬ 
tennas (Ef, Jbl, Wb and On) resulting in a maximum base¬ 
line of 7.5 MA (resolution of 27.6 mas). This allowed a reduc¬ 
tion in the number of channels to 32 per sub-band and in¬ 
creased time averaging to 5 seconds, while keeping the total 
jaflaplltude loss due to time and bandwidth smearing < 5%. 
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Moreover, the decrerised resolution lowered the number of 
pixels required to image a radius of 21" per phase centre 
and, therefore, the data size per phase centre, thus reducing 
the computation required to image the phase centres. Dirty 
maps of all phase centres were produced via the aips task 
IMAGR, using parameters: cellsize = 2.7 mas; image size = 
16384 X 16384 pixels; and natural weighting. Finally, 512 
pixels were removed from the image edges, to reduce alias¬ 
ing effects. This resulted in image sizes covering 41” x 41” 
per phase centre. The aips task imean was then used to ob¬ 
tain the RMS noise levels across the entire image, allowing 
comparisons with the peak pixel brightness. 

The distribution of the peak S/N for the 192 phase cen¬ 
tres is displayed in Figure |4] Six phase centres were found 
to have a peak S/N > 7 (j, with the remaining phase centres 
distributed between 5a and 6.15cr. While the six peak S/N 
are indeed real detections (and are discussed in Section [5]), 
it is possible that there are real detections with S/N < 7a. 


3.2 Source Identification 

3.S.1 Defining the deteetion threshold 

Morgan et al. (2013) addressed the problem of determining 
the likelihood that a bright pixel in a large VLBI image is 
a source (rather than a Gaussian noise spike). Consider a 
VLBI image of pixel size ripixeis x ii-pixeis- Assuming that 
the noise in the image is Gaussian distributed and is free 
from RFI or emission from real sources, the probability that 
a pixel exists with S/N greater than some threshold, Sm, 
can be shown to be 


P(S/N>5^) = l-erf (1) 

where erf(x) is the error function associated with the 
Gaussian distribution. Equation [T] is the probability that an 
image pixel of S/N is above a given threshold, assuming the 
noise is Gaussian distributed. Thus, pixels found above this 
threshold can be considered as being from real sources. The 
reciprocal of Equation[T]is therefore the equivalent frequency 
of occurrence, A/req, which is a measure of the total num¬ 
ber of image pixels such that a single pixel will have S/N 
equivalent to Sm- 

Figure [S] plots the pixel distribution of a single 15360 
pixels X 15360 pixels wide-held VLBI image, representing 
a blank region of the sky, 5 arcminutes from the pointing 
centre. The held covered by this im age was found to b e de- 
void of radio s o urces in the FIRST (iB^ker et al.|[l995h and 
iMaddox et al.l (l2007l l radio catalogue^, and lacked pixels 
with a hux density above the Gaussian noise threshold. The 
distribution follows a Gaussian distribution (the solid line), 
suggesting that the pixels of the wide-held VLBI images are 
Gaussian distributee^ 

Since multiple phase centres provide a method to image 
a large held through Vfleid smaller helds, the total number 
of pixels is Vfleid x ripixeis X ripixeis- By assuming the pixel 


® See Section IS.ll for a description of these radio source catalogues 
® Note that this is only valid if the data are free from real sources 
and RFI (a chievable through ca reful flagging). This is shown em¬ 
pirically in [Morgan et all (1201311 and Figure |5] 
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Figure 5. Histogram showing the pixel distribution of a single 
15360 pixels X 15360 pixels wide-field VLBI image, with the y axis 
in log scale. The solid line represents a least squares Gaussian fit 
to the distribution, with the mean and standard deviation given 
inset. The dashed lines plot the location of the 6.7 (t detection 
threshold. 



brightness distribution in the wide-held images is Gaussian 
distributed, from Equation [T] one 5.8 (t pixel would be ex¬ 
pected on average in every 15360 pixel x 15360 pixel image 
and one 6.7a pixel in 192 images, assuming no astronomical 
sources. This is in agreement with Figure [S] which clearly 
shows there are no pixels on either ends of the distribution 
beyond 6.7a. Thus, any peak hux density with S/N ^ 6.7a is 
more likely to be a source than a noise spike (assuming real 
sources are the only thing that could make the distribution 
non-Gaussian). 


3.2.2 Using Prior Catalogues to Search for Faint Sources 

The detection threshold calculated in Sect. 13.2.11 assumes 
that the image pixels are independent. This is not the case, 
as pixels within an area dehned by the synthesised beam may 
be correlated. With a cellsize of 2.7 mas and resolution of 
27.6 mas, gives ~100 pixels per synthesized beam, reducing 
the threshold to 5.1 ct in a single 15360 pixel x 15360 pixel 
image and one O.Oa pixel in 192 images. 

However, there is a possibility that real sources exist 
in the wide-held images below this limit and are treated as 
a noise-spike. To determine thi s probability the like lihood 
ratio (LR) method described bv [Morgan et al.) l|2013h is ap¬ 
plied. This method calculates the likelihood that each pixel 
in a VLBI image is associated with a catalogued source of 
known position, size and hux density (the alternative hy¬ 
pothesis, Hi), or simply a noise spike (the null hypothe¬ 
sis, Ho). This calculation is made assuming the noise in 
each image is Gaussian distributed and the RMS of the 
noise does not vary across an individual image. This method 
generates LR images corresponding to the input VLBI im¬ 
ages, where the pixel values correspond to the likelihood 
that the p ixel is related to a catalogued source (see Sec¬ 
tion 2.5 of [Morgan et al.l [20131 for a detailed description of 
this method). 

We use this method to search for VLBI counterparts 
of sources detected in M51 by the VLA FIR ST survey 
([Becker et al.l [l995l) and [Maddox et al.[ |[2007l) using the 
VLA in A configuration. In their analysis, Morgan et al.l 
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ll2013ll found that the results of the LR test did not give 
a significant advantage over simply detecting bright pixels 
due to the poor constraints on the position given by the 
lower resolution catalogues they used. Given t he better po¬ 
sition al accur acy provided by the VL A FIRST dBecker et al.l 
Il99flh and the lMaddox et al] |2003) surveys, it was deemed 
worthwhile to attempt to use this method. 

However, similar to [Morgan et al.l ll2013ll we found that 
the LR of a bright pixel was more dependent on the VLBI 
pixel S/N (i.e. p(H’o)), than the positional constraints given 
by the lower resolution catalogues (see Appendix A for fur¬ 
ther details on this method and the results). As such no 
significant advantage is given by the LR method over sim¬ 
ply detecting the brightest pixels. Thus, it is difficult to 
determine whether pixels with S/N< 6.7 (t are indeed real 
detections, and are henceforth considered as resulting from 
noise. 


were required (e.g uniform illumination, symmetric primary 
beam, and a perfect parabola defined by a Bessel function), 
which would no doubt affect the final correction. While there 
are errors in the final primary beam-corrected flux values, 
they are considerably less than not applying the primary 
beam correction. Thus, there is no prospect of discovering 
new sources by re-imaging with a primary beam correction 
applied to the visibilities. Nonetheless, in estimating the flux 
of the sources, the overall primary beam correction is taken 
into account and the source fluxes are adjusted accordingly. 

For radio sources within the primary beam of the Ef 
antenna (radius = 3.8') we estimate the residual amplitude 
error due to the primary beam is on the order of 10%. While 
for sources beyond the HPBW of the Ef primary beam, an 
error of 20% is estimated. This term is added in quadrature 
with the other sources of error which is listed in Table [T] 


3.4 Flux Density and Morphology 


3.3 Primary Beam Correction 


Primary beam corrections in VLBI observations are 
uncommon and to-date hav e been attempted for a 
subset of EVN antennas (|Mprgan__eLaL 2011 ) and 
the VLB A ( Middelberg et ed] 201 ll : Morgan et al. l2013l : 
iMiddelberg et al.ll2013l l. Unlike interferometers such as the 
JVLA where the primary beam correctio n is performed in 
the im age plane, the approach taken by IMiddelberg et aH 
l|201lh determines the primary beam attenuation through 
calculating the visibility gains. This requires the assump¬ 
tion that the correction is constant across the entire image. 

Primary beam corrections were applied to the phase 
ce ntres containing p i xels a bove 6.7 (t, following the m ethod 
of IMiddelberg et aP 1 I 2 OIII I and [Morgan et aLl 1 I 2 OI 3 I I . Very 
limited information is available on the primary beams of 
EVN dishes and so our primary beam model is necessarily 
simplistic, assuming a rotationally symmetric Bessel func¬ 
tion with a EWHM calculated from the published diameter 
of the disl0 for both polarisations. Eor Effelsberg only, we 
used the EWHM quoted in[3 scaled relative to the reference 
frequency quoted therein. Eor all dishes the primary beam 
was calculated using the central frequency of each subband, 
for the duration of the observation, using the measured sys¬ 
tem temperatures for each antenna and then applied to the 
visibilities. 

It should be noted that the primary beam correction 
did not improve the S/N of the detected sources. This is 
evidence of fairly signihcant errors in the model, since a cor¬ 
rect model would result in optimal weighting of each base¬ 
line. For the EVN, unlike the VLBA array, the difference in 
antenna sizes (from 100 m to 25 m) results in very different 
primary beam responses. Moreover, during the course of the 
observations the different physical characteristics of the indi¬ 
vidual elements produce varied effects on the overall primary 
beam response of the arrajQ- Furthermore, the exact char¬ 
acteristics of each antenna in the EVN (except for Ef to a 
certain degree) are unknowiF^ and a number of assumptions 


In total si x sources were fo und using the blobcat source 
extractor llHales et al.ll201^ ) with flux densities above 6.7(7, 
which are listed in Table [T] contour maps and locations of 
the sources in the field are displayed in Figures [6] and [71 
respectively. Following the correction for the amplitude re¬ 
sponse of the primary beam, the positions of the six sources 
were shifted to the image centre via UVFIX, averaged in fre¬ 
quency, and reimaged. 

CLEANed images were obtained using data from the full 
array via the aips task imagr. The images were made with 
cellsize = 0.5 mas, imsize = 4096 pixels x 4096 pixels, and 
100 clean iterations. Briggs robust parameter 4 was used to 
weight the images as it was found to produce an acceptable 
compromise between sensitivity and resolution. The final im¬ 
ages are displayed in Figure (6] 

Measurements of flux density and morphology were ob¬ 
tained by fitting a single 2D Gaussian component via IM- 
FIT, with independent measures of flux density from blob- 
cat and source sizes from fitting a cir cular Gaussian to the 
visibilities via modelfit in difmap (IShenherdl 1 1997ll . The 
measurements are list ed in TablelH with any c ross identihca- 
tion with the FIRST (iBecker et ^|[l995ll and iMaddox et H] 
(l2007tl surveys. The different measurements of flux density 
and source size were found to be consistent, so only the IM- 
FIT values of flux density and source size are included in 
Table [T] 

Furthermore, to search for possible emission associated 
with extended structures, lower resolution maps with the 
full array were made with a (u, v) range taper of 10 MA for 
all sources and are included in Figure |6] (the red contours). 
Each source will be discussed in detail in Section 5.5. 


4 MULTI-WAVELENGTH ANCILLARY DATA 


To fully understand the nature of the compact sources de¬ 
tected in the wide-held VLBI images, multi-wavelength data 

see http: //www. evlbi. org/user_guide/EVNstatus. txt 

https://efflOOmwiki.mpifr-bonn.mpg.de/doku.php?id=informatioii_for_astronomers;rx:p200imii 

Note:- this is also true for the VLBA beams of the EVN and MERLIN antennas, (R. Beswick, private 

There are current efforts being made to measure the primary communication) 





















































A High Resolution Wide-Field Radio Survey of M51 7 




(a) J132952+471142 


(b) J133005+471035 




Figur© 6. Sources detected with the full array of the EVN at 18 cm. The black contours are full resolution images with robust weighting 4, with a beam size 
= 9.86 mas X 5.92 mas and position angle = 14.3“^ . The lowest contours (for all but J133005+471010) are at it3<T, and increase in multiples of \/2, where a is 
the RMS of the images found via imean: 9.1 ;i.Jybeam“^ (J132952+471142); 17.4 //Jybeam“^(J133005+471035); 18.3 ;i.Jybeam“^ (J133011+471041); 
18.0 /.iJybeam”^ (J132932+471123); 18.6 ^Jybeam”^ (J133016+471024). The contours of J133005+471010 begins at zbl0% the peak flux density 
(3.26 mjy/beam), and increase in increments of 20% the peak flux. The sources were also imaged with a (t/.,u) range taper of 10 MA (red contours), with a 
resulting a beam size = 25.3 mas X 21.2 mas and position angle = -4.9° . The lowest contours are ib3cr, and increase in multiples of \/2 for; J132952+471142 
(<7 = 10.5 p,Jybeam~^), J132932+471123 (<t = 21.8 p,Jybeam~^) and J133016+471024 (<t = 24.4 p,Jybeam~^). While the remaining sources the contours begin 
at ±10% the peak flux density (J133005±471035, 506.6 /zjybeam”^; J133005±471010, 4.27 mJybeam~^; and J133011±471041, 552.6 /.tJybeam“^), and 
increase in increments of 20% the peak flux density . 








8 H. Rampadarath et al 


Table 1. Sources Detected in M51 with the EVN at 18 cm. Columns: (1) Source name; (2) & (3) Source position; (4) Angular distance 
to the centre of M51. At 8.4 Mpc, 1' ^ 2.44 kpc; (5) Signal to noise ratio of the source in the wide-field source detection maps; (6) Flux 
density measured from i mfit using a single 2-D Gaussian component; (7) & (8) Source sizes measured via imfit; (9) Cross identification 
with the lMaddox et ^ ll2007ll survey; (10) Detected by FIRST (Y), or not (N). 


Source 

(1) 

Position (J2000.0) 

0 

(') 

(4) 

S/N 

(5) 

Fimfit 

(^Jy) 

(6) 

IMFIT Size 

Identification 

Maddox FIRST 
(9) (10) 

RA 

13'* 

(2) 

Dec 

47° 

(3) 

Major 

(mas) 

(7) 

Minor 

(mas) 

(8) 

3132952+471142° 

29™ 52.7102(3)*’ 

11'42.7465(5)" 

0 

9.7 

140.1+16.0 

8.8+1.1 

1.9 +0.9 

53 

Y 

J133005+471035 

30™5.1381(1)*’ 

10' 35.7896(1)" 

2.39 

43 

2041.4+204 

12.8+1.9 

9.3+1.9 

104 

Y 

J133005+471010t 

30™5.1057(1)*’ 

10' 10.924(1)" 

2.60 

269 

4306.7+431 

6+0.1 

- 

- 

N 

J133011+471041 

30™11.0087(1)*’ 

10' 41.1393(2)" 

3.27 

32.7 

559.7+63.9 

6.1+3.0 

- 

107 

N 

J132932+471123 

29™32.1262(4)*’ 

11'23.1175(7)" 

3.52 

7.1 

188.1+51.5 

8.4+2.2 

3.1+1.5 

- 

N 

J133016+471024 

30™ 16.0265(1)*’ 

10' 24.3525(2)" 

4.17 

20.8 

455.2+64.5 

5.1+0.6 

- 

- 

Y 


*at M51a Centre; tSN 2011dh 



Figure 7. Location of the six sources detected with the EVN 
shown in Figure |6] and listed in Table [T] Flat spectrum sources 
(—0.5 ^ a ^5 0.5) are the pink circles and steep spectrum sources 
a ^ 0.5 are the green squares. The yellow elliptical source is 
the supernova SN 201 Idh, that exploded 159 days prior to our 
EVN observation (see Section 5.5.3). The green cross is a steep 
spectrum source detected with MERLIN at 18 cm (designated 
Component N) but not with the EVN at 18 cm (see Section 5.5.2). 

were obtained. Here, the data from different facilities that 
were used to complement the VLBI results are described, in¬ 
cluding any data processing and analysis that was required. 

4.1 Radio Continuum 

Archival radio data from the VLA and the Multi-Element 
Radio Linked Interferometer Network (MERLIN) were ob¬ 
tained for this study. M51 was observed by MERLIN on 
the 02 July 2005 at 18 cm. The data were calibrated fol¬ 
lowing the MERLIN user guidj^ The entire M51 system 
(11'.2 X 6'. 9) was imaged in aips using imagr, with Briggs 
robust parameter 4. Initial imaging of the full MERLIN ar¬ 
ray revealed a bright ~5 mjybeam”^ radio source coinci¬ 
dent with the VLBI source J133005-I-471035 and a weak de¬ 
tection above 300 /rjy beam“^ (i.e. 5a) at the centre of M51. 
By discarding the longest baselines (> 0.8MA) and contin¬ 
uum subtracting J133005-I-471035, deeper images were made 
revealing a number of radio sources, including diffuse, ex¬ 
tended emission around the core, previously detected only 

http://www.merlin.ac.uk/user_guide/OnlineMUG-ajh/ 


with the VLA. MERLIN sources coincident with the VLBI 
detected sources (Table [T]) were found via blobcat, using 
background RMS maps obtained from the source finding 
algorithm AEGEAN jHancock et al.l [20121 ). which were con¬ 
firmed via visual inspection. The final image gave an rms of 
50 /ijy beam~^and a beam size of 0.34" x 0.17" , PA=-6.8° . 

Archive VLA-B observations of M51 at 5 GHz (Obser¬ 
vation date of 11 December 2003) and VLA-A 1.4 GHz (Ob¬ 
servation date of 10 Oc tober 2004) were ob tained and cali¬ 
brated as described bv iDumas et al.l ll201ll ). Both datasets 
were imaged using the C ommon Astronomy So ftware Appli¬ 
cations (CASA) package dMcMullin et al.|l2007l) . with Briggs 
robust parameter 2, cell size of 0.3" and image size of 4096 
X 4096 pixels. The imaging procedure also took into consid¬ 
eration w-projection and finally corrected for the primary 
beam. The final image gave an rms of 12 pjybeam~^and 
a beam size of 1.08" x 1.02", PA=-87.2° for the 1.4 GHz 
observation and rms of 12 /ijybeam“^and a beam size of 
1.72" X 1.28" , PA=-79.6°for the 5 GHz observation. These 
results will be discussed in the individual source discussions 
in Section [S] 

4.2 Optical 

Optical images of the M51 system were acquired from the 
Hubble Space Telescope (HST) online archive. A six-pointing 
mosaic of the galaxy pair was obtained with the Advanced 
Camera for Surveys (ACS), Wide Field Channel (WFC) on 
board the HST in four filterQ B (F435W), V (F555W), 
I (F814W), and Ha R (F658N) lllJbedall2014l) . The ACS 
pipeline-reduced images (with bias subtraction, flat-fielding 
and drizzling applied) were downloaded from the MAST 
archive. The final surface brightness of the images are given 
in units of electrons per second, es“^, which can be converted 
to flux density by the expression 

• 3631 - 10 "''^''®""‘^'/^®', ( 2 ) 

where, is the flux density (in Jansky) at the ef¬ 

fective optical wavelength, Zab is the zeropoint and de¬ 
fines the absolute physical flux density for a point source 

HST observing program 10452, P.I. Steven V. W. Beckwith 
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Table 2. Parameters of the HST ACS Filters 


Table 3. Log of the chandra observations 


Filter 

Colour 

-^eft 

(A) 

^AB 

(mag) 

(mag) 

F435W 

B 

4450 

25.67 

0.125 

F555W 

V 

5510 

25.72 

0.091 

F814W 

I 

8060 

25.94 

0.053 

F658N 

Ha 

6563 

22.77 

0.08 


*^Siriann^_e^_^^005; ^ Obtained from the 
NASA/IPAC Extragalactic Database 


of 1 es~^ for a particular ACS-WFC filter llSirianni et alj 
l2005ll . and A\ is the Galactic foreground colour magnitude 
extinction towards M51. The values of Zab and A\ for the 
different HST filters are listed in Table [5] The parameter 
Cgs-i is the net count rate for an infinite aperture. We used 
standard Heasoft tasks to carry out aperture photometry 


by 

by 

Fx 


Sirianni et al. 

(20051 

Sirianni et al. 

to 

o 

o 

m 


and the Johnson-Cousins UBVRI photometric sys¬ 
tems. 


4.3 X-rays 

M51 has been observed by Chandra/ACIS-S 12 times (con¬ 
sidering only observations longer than 2ks) between 2000 
and 2012 (Table [3]). Aim points and roll orientations dif¬ 
fered between observations, but the nucleus of M51a was al¬ 
ways located in the field of the back-illuminated S3 chip. We 
downloaded the data from the public Chandra archive. We 
reprocessed each observation usin g tools available in CIAO 
Version 4.5 llFruscione et al.ll2006l l , and CalDB version 4.5.3, 
to obtain new level-2 event files. We then reprojected the 
event files to a common tangent point using the CIAO script 
reproject_ all. The total exposure time in the merged file 
is 855 ks for most of the D25 region of M51a. We used the 
merged event file for imaging and source detection purposes. 
For spectral analysis of the brightest sources, we generated 
separate PHA, RMF, and ARF files for each observation 
with the CIAO task specextract, and created combined spec¬ 
tra by averaging and merging the spectra from the individ¬ 
ual observations. We then fitted the spectra w ith standard 
models in XSPEC Version 12.6 (lArnaudlll996l l. For sources 
with sufficiently high S/N, we binned the spectra to 15 
counts per channel and used statist ics. For sou rces with 
< 100 counts, we used Cash statistics llCas 311^). 

Table |4] lists the X-ray detections for the sources in Ta¬ 
ble [T] and other sources discussed in this paper: SN 1945A; 
SN 19941; SN 2005cs; and Component N. 


ObsID 

(1) 

Obs. Date 
(dd/mm/yy) 
(2) 

Exp. Time 
(ks) 

(3) 

354 

20/06/00 

14.86 

1622 

23/06/01 

26.81 

3932 

07/08/03 

47.97 

12562 

12/06/11 

9.63 

12668 

03/07/11 

9.99 

13812 

09/09/12 

179.2 

13813 

12/09/12 

157.46 

15496 

19/09/12 

40.97 

13814 

20/09/12 

189.85 

13815 

23/09/12 

67.18 

13816 

26/09/12 

73.1 

15553 

10/10/12 

37.57 


Table 4. Radio/X-ray 

ASSOCIATIONS 




Position (J2000.0) 



RA 

Dec 


VLBI Sources 

lO'^ 

47° 

X-rays? 

J132952+471142* 

29'"52®.7102 

11' 42" .7465 

YES 

J133005+471035 

30'"05®.1381 

10' 35” .7896 

YES* 

1133005+471010* 

30’"05M057 

10' 10" .924 

YES 

J133011+471041 

30"* 11®.0087 

10'41” .1393 

YES 

J132932+471123 

29"*32®.1262 

11'23”.1175 

YES* 

J133016+471024 

30"* 16® .0265 

10' 24" .3525 

YES 

Other Sources 

RA 

Dec 

X-rays? 

SN 1945A 

29"* 58® .95 

15' 53" .8 

NO 

SN 19941 

29™ 54® .072 

11' 30" .5 

YES 

SN 2005CS 

29™ 52® .850 

10' 36" .3 

unclear 

Component N 

29™ 51®.59 

12' 07" .7 

YES* 


^Nuclear source; 1-SN 2011dh; * very faint 


5 DISCUSSION 

5.1 Fraction of Radio Sources Detected 

Of the 107 radio sources detected in M51 by the 
iMaddox et al.l ll2007l l VLA survey (beam size: 1.50” xl.2l” , 
21 cm; 1.47” xl.13”, 6 cm. Sensitivity: 22.5/iJybeam“^, 
21 cm; 11.7/iJybearn”^, 6 cm), we report three confirmed 
detections. Similarly, of the 15 sources in the same regio n 
detected by the VLA FIRST survey llBecker et al.|[l995li . 
three were detectecf^. For a more realistic estimate of the 
proportion of detected sources we restrict ourselves to the 
catalogued sources that are bright enough to be detected 
above the 6.7cr threshold of the EVN observation. The cat¬ 
alogued flux densities were adjusted to 18 cm by using the 
spectral indices provided by the catalogues, and the sensi¬ 
tivity of the EVN phase centres were correct ed for primary 
beam effects. This leaves 70 sources from the lMaddox et aP 
ll2007ll VLA survey and 15 from the FIRST survey. Assum¬ 
ing no catalogued sources are resolved, the fraction of VLBI 


Two sources were detected by both surveys - the core of M51 
and source J133005+471035 
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sources detected compared to the lMaddox et all ll2007ll and 
VLA FIRST surveys are 3/70 and 3/15 (tom) Q 

respectively. 

Considering sources from the lMaddox et aT] ll2007l l sur¬ 
vey not associated with HIT regions, the detection rate is 
3/32 (jlo'o 3 )i while the VLBI detection ra te for sources 
with peak flux density > 1 mjy in the iMaddox et al.l 
1 I 2 OO 7 II survey is 2/5 (Iq i 2 )- The fraction of sources de- 
tected from the FIRST survey (3/15) and iMaddox et al.l 
(l2007tl sources with flux density > 1 mJy (2/5), closely 
agrees with the d etection rates of prev i ous wide-held 
VLBI surveys (e.g. iM iddelbe rg et akl I 2 OIII : iMorgan et al.l 
201^^ Middelherg_etaL 20131) . However, the fraction of 


Maddox et al.l (l2007l l sources not associated with HII re¬ 


gions detected with VLBI (3/32) is significantly lower than 
the previous wide-held VLBI studies. 

At 1.6 GHz, our observations would be less a f fected 
by free-free absorption co mpared t o IM addox et al.l (l2007l l 

1.4 GHz observations fe.g iTingavl l2004h . It is likely that 
the majority of th e radio sources in M51 detected by 
IMaddox et al.l ll2007l l are fully resolved at our VLBI resolu¬ 
tion, suggesting perhaps either lack of compact radio emis¬ 
sion from sources associated with stellar clusters or the ex¬ 
pansion of supernova remnants (SNRs) beyond the maxi- 
mu m size detectab l e with our observations. 

IMaddox et al.l 1 I 2 OO 7 I ) estimated that the typical sizes 
for SNRs in M51 range between 0.3” and 0.5”, (i.e. 12.2 - 

20.4 pc). Gonsidering the maximum source size resolvable 
with our observations is 0.14” (~6 pc; see Section 12.111 . the 
non-detections are not surprising. Assuming an expansion 
veloc ity of 21,000 km s ~^ for core-collapse s upernovae in 
M51 dWeiler et al.ll20jT] : lBietenholz et al.ll201 j ) gives an age 
of 570 - 950 yr for the largest and 280 years for the smallest 
shell sources. 


5.2 The M51a Nuclear Region 

5.2.1 A Faint Seyfert Radio Nucleus 

The VLBI observation resolves the radio source at 
the nucleus of M51a into multiple components (source 
J132952-I-471142 in Figure |S^). The main feature is a com¬ 
pact source (labelled Gl) marginally resolved in the high 
resolution images (black contours), and unresolved in the 
10 MA image. Located to the north of Gl are two weaker 
sources (labelled G2a and G2b), which are blended into a 
single, unresolved elongated structure that appears to be 
connected to Gl in the 10 MA image. Gomponents G2a and 
G2b are unresolved sources with total flux densities of 51.5 
± 9.1 /rjy and 54.6 ± 10.2 /rjy, respectively. The separation 
of G2a and G2b from Gl are ~ 41.5 mas (~ 1.69 pc) and ~ 

49.4 mas (~ 2.01 pc), respectively. 

VLBI observations of a distance-limited sample 
(<22 Mpc) of Seyfert galaxies, at 18 and 6 cm, 
have re vealed similar stru c tures to M51a within their 
nuclei (iGiroletti &: Pa nepal l2009l : iBontempi et al.l I 2 OI 2 I : 
IPanessa fc Girolettil 120131) . Here, the detected VLBI flux 
density is dominated by a high brightness temperature 


Population uncertainties wer e derived using the Bayesian beta 
distribution quantile technique llCamero ^ l201l|) . 


(logioTfl > 7.5) parsec-scale source with a flat/intermediate 
spectrum (0.3 a 0.6), accompanied by fainter, steep 
spectrum extended components. The flat/intermediate spec¬ 
trum, high Tb sources are taken as evidence for non- 
thermal processes driven by jet-producing central engines 
typical of radio-quiet and radio loud qu asars and AGNs 
dBlundell fc Beaslevlll99^ : lDoi et al.ll20lij) . 

Using the brightness temperature relation, Tb 

— - f~') I where: S{v) is the flux density at 

z/Ct/majt/inin ^ U ^ 

frequency w, k - the Boltzmann constant; ^maj and^min are 
the major and minor deconvolved sizes in radians; and c 
is the speed of light, a value of logioTs ~ 6.6 is obtained 
for Gl. Both Tb and luminosity (lO*^® W Hz“^) of Gl 
are within the typical range o f values for low luminosit y 
AGNs within Seyfert galaxies dPanessa fc Girolettil |201^ . 
Therefore, components G2a and G2b are, most probably, 
hotspots associated with a jet or thermal outflow. 


5.2.2 Previous VLBI observations 

iBontempi et al.l 1 I 2 OI 2 I ) observed M51a with the EVN at 
both 18 and 6 cm and were unable to detect any emission 
above 3a (75 and 160 /ijybeam“^ at 18 and 6 cm, respec¬ 
tively). However, they searched within a 1” radius centred 
on the VLA FIRST position RA = 13^ 29™ 52.804°, Dec = 
4-47° 11 ' 40.065” , which would mean that the VLBI position 
in Table [T] was beyond their search field. The uncalibrated 
datasets for both observations (Experiment ID: EG037D 
and EG037G. PI: Giroletti) were obtained from the EVN 
data archive. Using the calibration tables generated by the 
EVN pipeline, both datasets were imaged as described in 
Section 5.3.1. 

The core-jet structure was detected only in the 18 cm 
data at the expected position with a similar size, flux den¬ 
sity and morphology as seen in the wide-held observation 
presented in Section [3.41 The non-detection at 5 GHz is at¬ 
tributed to the poor sensitivity of the observation (ct= 50 
/rJybeam”^), caused by the loss of Ef for 90% of the 
time due to high winds. Assuming a hat to intermedi¬ 
ate spectrum (—0.5 ^ a ^ 0.5), the expected hux den¬ 
sity of J132952-I-471142 at 6 cm would be less than 150 
pjy beam~^. 

The images from both epochs (2008 and 2011) at 18 cm 
are compared in Eigure [S] where the contours of the 2011 
epoch are overlayed on the 2008 epoch (colour image). A 
restoring beam of size 18.5 mas x 11.0 mas with a P.A. = 
-3.71° was used for the 2011 epoch so as to match the 2008 
epoch. While the compact radio core is found at the same 
position in both epochs, the lower brightness compact source 
associated with a jet/thermal outflow shows an apparent 
lateral shift in position of 7.53 mas between epochs. The 
positional measurementsfor both epochs were obtained using 
the same phase reference source, with the same referenced 
position. 

With a time difference of 3.7 years, and assuming that 
the difference in position is genuine, a velocity of 0.27c is im¬ 
plied. Eurthermore, the peak brightness of the jet/hotspot 
has decreased between epochs (80.2 /rjybeam“^, 2008; 44.2 
pjybeam~^, 2011), which may be associated with the ap¬ 
parent motion. Molecular aperture s ynthesis HGN observa¬ 
tions of M51a bv iKohno et aP lll996l) revealed a twisted ve- 
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13*’29“52.714® 52.712® 52.710® 52.708® 52.706® 52.704® 

Right Ascension (J2000) 


Figure 8. VLBI images of t he Seyfert nucleus of M51a for the 
2008 epoch (colour) with the lBontempi et al.l ll2012li dataset and 
the 2011 epoch (contours). Both epochs were imaged with a beam 
size of 18.5 mas X 11.0 mas with a P.A. = -3.71° . The RMS noise, 
a in both images are 19.2 ^Jy beam“^ and 8.9 /rjy beam“^ for the 
2008 and 2011 images respectively. While the peak fluxes are 102.3 
it 21.0 ^Jybeam~^ (2008) and 84.7 it 14.7 ^Jybeam“^ (2011). 
The first contours of the 2011 image are at it3cr, and increase in 
multiples of \/2. The offset in position of the low surface bright¬ 
ness structures between epochs is 7.53 mas. A third EVN obser¬ 
vation has been requested to investigate the nature of this offset. 


locity field resulting from the existence of a dense rotating 
molecular disk that is aligned to the nuclear jet of M51a, 
and the circular rotation of the galactic disk plane. How¬ 
ever, the maximum velocity of this field (450-500 km s“^) is 
too low to result in the apparent motion of the pc-scale jet. 
Furthermore, the direction of the apparent motion suggests 
that the jet may be precessing. While j et precession is not 
uncommon in radio bright AGNs (e.g IConwav fc MurphvI 
I 1993 I : IAppI et al.lll99d : ILu fc Zhoull2005l ). it is vet to be ob- 
served for low luminosity Seyfert galaxies. In addition, out¬ 
ward proper motions have been observed (at the 3cr level) 
in the jets of the Seyf ert galaxy NGC 5506 on pc scales 
llMiddelberg et al]|2004h . The apparent motion could be the 
result of a combined effects of forward bulk m otions of the 
j et int eracting with the rotation observed by iKohno et al.l 
lll996li . 

It is, however, possible that the position shift can be 
the result of systematic factors (e.g. poor calibration), since 
the hotspots displaying the apparent motions are ~4 (t de¬ 
tections. This is unlikely since the positions of the phase 
calibrator and Cl between both epochs were found to agree 
to within less than 1 mas. Whether this apparent motion is 
indeed intrinsic to M51a is difficult to ascertain with only 
two epochs. Therefore, a request has been made to observe 
the nuclear region of M51a with a third higher sensitivity 
observation with the EVN at 18 cm to determine whether 
this is indeed intrinsic to the parsec-scale jet of M51a, the 
results of which will be discussed in a later paper. 


5.2.3 A Multi-wavelength Perspective 

Figure [9] shows the complex nuclear region of M51a with 
increasing resolution. The VLA-A image at 20 cm (reso¬ 


lution ~1.08" ) and the higher resolution MERLIN 18 cm 
images (resolution ~0.34" ) reveal an area of ring-like emis¬ 
sion to the north of the core (hereafter, the northern bub¬ 
ble) and a dense area of emission to the south, referred 
to as an extranuclear cloud (XNC), that may have re¬ 
sulted fr om interactions of a nuclear jet with the surround¬ 
ing IS M llCrane fc van der Hulst||l993;l Terashima fc WilsonI 


120011 : iMaddox et al. 2007 ). Maddox et al.l i 20071) mea¬ 
sured peak flux densities for the nucleus at both 20 cm 
(2.08 mjybeam“^) and 6 cm (1.14 mjybeam“^), giving 
a moderately fiat spectrum (a = —0.49). The MERLIN 
image shows the more compact regions where the ejecta 
are strongly interacting with the surrounding medium. Also 
shown in the image is a slightly extended source located 
~30” (~1.22 kpc) from the nuclear radio source, which is 
discussed later in Section (5.2.41 


As noted by previou s studies dTerashima fc WilsonI 
I 2 OOII : iMaddox et al.l 120071 ). the radio continuum emission 
of the nuclear region, including the XNC and the northern 
bubble is very similar to t he Chandra X-ray morphology. 
iTerashima fc WilsonI (1200 ll ) found that the X-ray spectrum 
mainly consist of a hard component that dominates the 3- 
8 keV band corresponding to an obscuration of the nuclear 
source by a column density in excess of 10^"^ cm~^. With 
our deeper combined spectrum (855 ks), we confirm that 
the emission consists of a soft thermal-plasma component 
plus a hard component probably due to reflection or scat¬ 
tering, with a strong Ee line. The average observed flux in 
the 0.3-10 keV band is fo.3-10 = (1-4 ± 0.1) x lO”*^^ erg 
cm~^ s“^, of which « 3.5 x lO”*^"^ erg cm“^ s“^ is in the 
thermal-plasma component, which may come from a com¬ 
pact nuclear starburst or from gas shock-ionized by nuclear 
outflows. The unabsorbed luminosity of the thermal-plasma 
emission is Lo.3-10 ~ 3.3 x 10^® erg s“^. For the hard com¬ 
ponent, nominally we obtain Lo.3-10 = (9±1) x 10^® erg s“^; 
however, this is only a minimal fraction of the true X-ray 
luminosity of the nuclear black hole, whose direct emission 
is complete ly obscured in the Cha ndra energy band. Using 
BeppoSAX, iFukazawa et al.l ll200ll ) detected the direct hard 
X-ray emission at energies > 10 keV, and estimated a col¬ 
umn density Ah « 6 x 10^"^ cm“^ and an intrinsic nuclear 
luminosity L2-10 ~ 10"*^ erg s“^. More detailed analysis of 
the nuclear X-ray spectrum is beyond the scope of this work. 

The X-ray emission of both the XNC and the north¬ 
ern bubble is dominated by optically-thin thermal plasma 
at a characteristic temp erature kT ~ 0.6 keV, ty pical of 
collisionally-ionize d gas ( Terashrmn^_j£_Wilson|l200ll ) . High- 
resolution 6 cm llCrane &: van der Huls^ 19921) and 3 cm 
dBradlev et al.|[2004h VLA observations of the nucleus re¬ 
veals a narrow jet-like feature emanating from the nucleus 
and terminating in the XNC. [O ill] line measurements of 
narrow line regions within the nuclear region, and along the 
direction of the narrow j et-like feature , find evidence for 
shock heating of the gas (iBradlev et al.l [20041 ) at the loca¬ 
tion of the XNC, indicating continuous fuelling of the XNC 
by the radio jet. However, no evidence has been found for 
shock ionisation in narr ow line regions with in the direction 
of the northern bubble dBradlev et al.ll2004h . suggesting the 
northern bubble may therefore be f rom a previous ejection 
cycle of the Seyfert nucleu s of M51a dCrane fc van der HulstI 
I 1992 I: IMaddox et al.ll200^ . Interestingly, jet hotspots were 
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VLA-A 21 cm MERLIN 18 cm 



Right Ascension (J2000) Ascension (J2000) 


Figure 9. Radio images of the nuclear region, XNC, northern bubble and a bright source located ~ 30” (~1.22 kpc) from the core. 
VLA-A: The first contours are at ±15% the peak flux (11.49mjy beam“^, and increase in increments of 10% the peak flux; MERLIN: 
The first contours are at ±3 (t, and increases in multiples of \/2, where cr = 50fj,Jy beam“^; EVN: see Figure [fil.. The negative contours 
are plotted in grey. 


detected with VLBI only in the direction of the northern 
bubble (components C2a and C2b). However, the AGN core 
appears slightly elongated in the VLBI image, in the direc¬ 
tion of the southern jet, and may be the launching point for 
the jet. 


5 . 2.4 Component N: a possible fossil radio hotspot 

The VLA and MERLIN images reveal a bright object 
(2229 ±21 pjybeam“^, VLA 20 cm; 946 ±11 ^Jybeam“^, 
VLA 6 cm; 448.3 ± 54 /ijybeam“^, MERLIN 18 cm;) lo¬ 


cated to the north of the core of M51a that was not de¬ 
tected in the EVN observations. The MERLIN position (RA 
= IS** 29“ 51.593°, Dec = 47° 12' 07.66” ) places Component 
N ~25.85” from the nearest phase centre. The unaveraged 
visibilities for the phase centre closest, was shifted to the 
position of Component N. The dataset was then averaged in 
time and frequency, such that the combined smearing at a 
radius of 2” from the MERLIN position was limited to 5%. 
The resulting dataset was imaged also with uv-tapering of 10 
MA, with natural weighting, which resulted in no detection 
above 5a. Given that the primary beam sensitivity at com- 
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ponent N is is 5-10% the peak sensitivity, we estimate an up¬ 
per limit of the EVN peak flux density of 82.5 ^Jybeam“^. 

The object, hereafter component N, has a steep spec¬ 
trum (aicm" = "0.7), dominated by synchrotron radio emis¬ 
sion, thus suggesting it is possibly a jet hotspot associated 
with the Seyfert nucleus of M51a, a SNR or a background 
AGN. The MERLIN 18 cm image in Eigure|9] reveals a com¬ 
plicated structure, elongated in the direction of M51a nu¬ 
cleus, atypical of SNRs. Moreover, the lack of detection of 
high brightness temperature emission with the EVN sug¬ 
gests it may be a background AGN whose core has either 
been scatter broadened or a background galaxy without a 
compact radio core. 

There is a faint, point-like X-ray source IFigure fTOl top 
right panel) at the position of Component N in the Chandra 
data. The source is significantly detected in the hard band 
(2-7 keV). The detection is only at the 3a level; however, 
the exact positional coincidence with the radio source gives 
us confidence that it is a real X-ray source. At softer ener¬ 
gies, such a faint source is swamped by the diffuse thermal- 
plasma emission in that region (Figure [TITl top left panel); 
but we cannot rule out the possibility that the point-like 
X-ray source is heavily absorbed below 2 keV, perhaps by 
one of the numerous dust lanes in that region. HST contin¬ 
uum images do in fact suggest the presence of a dust feature 
at the location of Component N (Figure 1101 bottom right 
panel). Based on the Chandra/ACIS 2-7 keV count rate of 
(2.0 ± 0.6) X 10“® ct s“^, and assuming a power-law photon 
index F = 1.7 and line-of-sight absorption, we estimate an 
observed 0.3-10 keV flux / 0 . 3-10 = I-IIq's x 10“^® erg cm~^ 
s“^ and an emitted luminosity Lo. 3-10 = I.OIq'^ x 10®^ erg 
s“^ if located in M51. 

We inspected the HST/ACS-WFC images but did not 
find any obvious counterpart of Component N in the con¬ 
tinuum (Figure fTOl bottom right panel); a few nearby point¬ 
like sources are consistent with small, unresolved (FWHM 
< 0 ” 1 « 40 pc) groups of young massive stars, abundantly 
scattered in the inner disk of M51a. An elongated stream 
of Ha emission near Component N is oriented towards the 
nucleus (Figure [TOl bottom left panel), roughly in the same 
direction as the putative jet, but we cannot determine at 
this stage if there is a physical relationship. Further analysis 
of the circumnuclear Ha and continuum emission is beyond 
the scope of this paper. 

Possible interpretations for the Component N radio and 
X-ray sources are synchrotron emission from a hot spot en¬ 
ergized by a nuclear jet, or a background radio galaxy seen 
through the disk of M51a. If it is a background AGN, it is a 
remarkable coincidence that its radio structure is elongated 
towards the core of M51a, but that is not a priori impossible 
(a similar situation is well known in the spiral galaxy M 83, 
with a bright background FRH radio galaxy projected near 
the nucleus, which might a t first sight be interpreted as a 
jet fe ature from M 83 itself (iMaddox et al.|[2006l : iLong et al.l 
l2014h . 


To determine the probability of a chance align¬ 
ment of a background AGN with the M51a nu¬ 
cleus we used the method o f estim ating quasar align- 
ment by Edmunds fc George! (Il985l l. In this method, 
lEdmunds &: George! lll985ll investigated the expected num¬ 
ber of random alignments, /r, in a field of area A, with a 


surface density of sources, n for random clustered and un¬ 
clustered fields and found it to be approximately, 

U ~ ? ( 3 ) 

where dmax is the distance between the two furthest 
sources, Pmax is the distance between the outer two sources, 
such that the product dmax • 2pmax forms a rectangle enclos¬ 
ing all points of interest. The number density of quasars 
from the FIRST s urvey is ~90 quasars per square degree 
dWhite et al.lll99^ and the total area surveyed in this study 
is 11.2' x6.9'. Thus, using Equation [3] fj, 1 for the in¬ 
ner region of M51a i.e. dmax = 30" and Pmax = 6" . There 
is, therefore, a very low probability of a chance alignment 
within the central region of M51a, suggesting that compo¬ 
nent N and the nucleus of M51a are most probably physi¬ 
cally related. Thus, component N is most likely a jet hotspot 
associated with the Seyfert nucleus of M51a. 

However, with no evidence of a continuous jet in the 
northern region of M51a, and the lack of a high brightness 
temperature component. Component N may be a fossil ra¬ 
dio lobe/hotspot associated with past AGN activity from 
the core of M51a. Fossil radio lobes ha ve been detected 
in a few radio galaxies such as 3C 388_l [jones_&^restonl 
I 2001 I : iGentile et al.1 120071 ') and IC 2476 llCordevlll987ll . with 
a fe w “dying galaxies" with evidence for f ossil radio lobes 
(e.g. IParma et ani2007l : TMurgia et al.ll20i3l . Radio lobes or 
hotspots are powered by energy from an AGN via jets of 
plasma. When the AGN undergoes a fading or dying stage, 
the jets and eventually the hotspots will expand adiabati- 
cally into the surrounding medium, and disapp ear due to the 
loss o f the continuous production of plasma (iMurgia et al.l 
I 2 OI 2 II . It is possible that this eventually leads to the loss of 
the high brightness temperature component, leaving behind 
a diffuse shell. 

If this is the case then the remaining synchrotron ra¬ 
diation in component N would be from electrons older 
than those within the core (e.g. ICarilli et alJ ll9^L Syn¬ 
chrotron ageing in radio galaxies is generally associated with 
a spectr al “break" that shif ts to lower frequencies with older 
plasma (ICarilli et al.|[l99lh . To investigate this possibility, 
the inner region on M51a (including Component N) will be 
observed with the Low Frequency Array (LOFAR) including 
all the international stations in VLBI mode at 110 and 170 
MHz. 


5.3 Supernovae 

At the time of our observations, M51a and M51b 
hosted four opti cally observed supernova e (SNe): the Type 
la SN 1945A_j[KrwaI_fc Sax geri^ Il97lll : the Type Ib/c 


SN 1994 1 jPuckett et al! 1994li : the Type Ib/c SN 2005cs 


I Muendlein et_alj 120051 ') : and the Type Hb SN 201 Idh 


I Griga et al.l 20l3). Of the four SNe, ra dio emission and X- 
rays have been detected in SN 19941 (IWeiler et al.l 1201 ll . 
Radio), ([immler et al.l 200j^ 19^ . X-R ays) a nd SN 201 Idh 
llKrauss et al.l 2012|: [Bietenholz et al.l |2012|: [Horesh et al. 
2013, Radio), Isoder berg et al.l I 2 OI 2 I : ICampana fc Immler 
2012 ; iMaeda et al.ll2014l . X-Rots'). 

SN 2005cs is located within a region with a number 
of faint sources near the detection limit. There is however. 
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Figure 10. Top left: Chandra/ACIS-S X-ray image of the nuclear region of M51a. Colours are; red = 0.3—1 keV, green = 1—2 keV; blue 
= 2—7 keV. The location of the radio nucleus (coincident with the X-ray nuclear source) and of MERLIN’S Component N are circled. The 
radius of the circles is 1” ; this is chosen for display purposes only, and is much larger than the positional uncertainty. Note the similar 
morphology of the diffuse soft X-ray emission with the diffuse radio emission (Figure 11). Top right: Chandra/AClS-S X-ray image of 
the same field, in the hard band only (2—7 keV). The image has been smoothed with a 2-pixel Gaussian. Nucleus and Component N are 
circled. Bottom left: HST/ACS continuum-subtracted image in the F658N filter (Ho) of the same field, with nucleus and Component 
N circled. Bottom right: HST/ACS image in the F555W filter (V band), with nucleus and Component N circled. 


soft X-ray emission inside a 1” error circle, from at least 
two sources with «20-30 counts, but it is unclear whether 
either of them is the true X-ray counterpart of SN 2005cs. 
In any case, we can say that the X-ray luminosity Lo, 3 _io < 
2 X 10^® erg s“^, which is two orders of magnitude lower 
than the upper lim its on the X-ray luminosity estimated by 
iBrown et alJ ll2007l ) . We will now discuss the radio and x-ray 
properties of, SN 19941 and SN 2011dh. 


5.3.1 SN 19941 

Radio Almost a decade after the explosion, iMaddox et al.l 
II 2 OO 7 I) detected radio emission from SN 19941. How¬ 
ever, no radio emission above the detection threshold of 
50.0 /ijybeam“^ was found in the VLBI images at the 
position of SN 19941. G iven that our observ ation was al¬ 
most a decade from the iMaddox et al.l ll2007l ) observation, 
it is possible that SN 19941 has faded below our sensitivity 
and/or expanded beyond the EVN’s maxi mum resolvable 
size. A djusting the 1.4 GHz flux densit y from IMa ddox et ^ 
ll2007l) with a spectral index of -1.04 ( Maddox et al. 20071) . 
the fl ux density of SN 19941 at 1.65 GHz at the time 
of the IMaddox et al.l ll2007l ) observation would have been 
134.9/rJy beam~^. The supernova light curve model of 
IWeiler et al.l 1 I 2 OO 2 I ) indicates that the flux density, S de¬ 
creases with time following a power-law, S oc where 

tagu is supernova age and /3=-1.42 for SN 19941 (IWeiler et al.l 


At the time of our observation (age of 6428 days) a 
1.65 GHz peak flux density of 59.5 /rjy beam“^ is estimated, 
which is just above our detection threshold. Furthermore, 
the expansion of SN 19941 has b een well studied and follows 
the form 1.29(tage/l day) /ras l| Weiler et al.l l201lh . which 
estimates a size of 8.3 mas at the time of our observation. 
The estimates of size and flux density are just within the 
limits of our EVN observation; as these estimates are model 
dependent, it is possible that SN 19941 was just below the 
detection threshold. It is also possible that the power-law 
decrease of the flux density steepens at later times, reduc¬ 
ing the flux density below our detection threshold. This has 
been observe d for SN1957D, al beit more than 40 years after 
its explosion llLong et al. I I 2 OI 2 I) . 

X-Rays. SN 19941 was detected in the first three Chan¬ 
dra observations (Table [ 3 ]) with declining flux: /0.3-10 = 
(2.0 ± 1.0) X 10”^® erg s^^ cm"^ on 2000 June 20 (t = 2271 
d after the explosion), fo.3-10 = (1-7 ± 0.7) x 10“^® erg 
s“^ cm“^ on 2001 June 23 {t = 2639 d), and fo.3-10 = 
(0.9±0.4) X 10“^® erg s~^ cm“^ on 2003 August 7 (t = 3416 
d). The spectrum is soft (/0.3-2 = (0.92 ± 0.02) x /0.3-10), 
and well fitted with optically-thin thermal plasma emission 
at kT = (0.9 ± 0.1) keV. The corresponding unabsorbed lu¬ 
minosity on 2000 June 20 was fo.3-10 = (1.8±0.9) x 10®^ erg 
s~®. Subsequent observations were either too short for reli¬ 
able measurements (ObsIDs 12562 and 12568), or the source 
had become too faint to be significantly detected in individ- 
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ual observations. However, after stacking the seven Chan¬ 
dra observations from 2012 September-October {t « 6750 
d), we detected it at an average observed flux fo. 3-10 = 
(4 ± 1) X 10“^® erg s“^ cm“^, corresponding to an emitted 
luminosity Ri 4 x 10®® erg s“®. 

The X-ray luminosity evolution at late times is a func¬ 
tion of the progenitor wind properties as well as the SN 
shell expansion velocity dFransson et al.l[l996l : llmmler et al.l 

I2OO2II : 



where [fimp) is the mean mass per particle (~ 2.1 x 10“^'* 
g for a H-fHe plasma), A{T) is the cooling function (« 
3 X 10“^® erg cm® at 1 keV), M is the mass loss rate 
in the stellar wind in the last few 10“^ years before the ex¬ 
plosion, Hw is the stellar wind speed (~ 10 km s“®), Vb is 
the shell expansion speed (for SN 19941, Vs « 16, 500 km 
s“®: Filippenko et al. 1995), and t the time after the explo¬ 
sion. In the simplest scenario, for constant A(T), (M/uw) 
and Vb, the X-ray luminosity evolves as Using three 
ROSA T observations, and the first two Chandra observa¬ 
tions, (llmmler et al.l 1^0021 . Fig. 2) found an X-ray flux de¬ 
cline /x oc t~‘‘ with 1 < s < 1.5, suggesting a mass-loss rate 
M ~ 10“® Mq yr“®. We can now better constrain the late¬ 
time decline rate, by adding the flux datapoints from the 
later Chandra observations. We find that /x oc which 

is similar to the power-law decrease of the radio flux density. 


5.3.2 SN 2011dh 


Radio. Our VLBI observations were conducted 159 days 
after the explosion date (2011 May 31, hereafter to) of 
SN 2011dh. The observations detect a compact unresolved 
source with a flux density of 4.306 ± 0.024 mjy at the po¬ 
sition of SN 201 Idh (source J133005-I-471010). 

The position of our peak fl ux density is offse t from 
the VLBI positions reported by iMarti- Vidal et aP 1I2OIII I 
at 22 GHz 14 days after to and iBietenholz et al.l 1 20121 1 at 
8.4 GHz 179 days after to, by 1.9 mas and 2.3 mas, respec¬ 
tively. The positional measurements were obtained using the 
same phase reference source, with the same referenced posi¬ 
tion. The offset in position with the higher frequency VLBI 
is less than half the synthesised beam of our 18 cm EVN 
observation. 

With EVN observations at 8.4 GHz 179 days after to, 
IBietenholz et al] 1 I 2 OI 2 I I estimated the radius of the shock 
front of SN 2011dh to be 0.25±0.08 mas. The radius of 
the shock front was also estimated up to 92 days after to 
via the radio spectral energy distributi on, fitted to a syn¬ 
chrotron self-absorpti on (SSA) model llKrauss et al.l l2012l : 
IBietenholz et al]l2012ll . Combining the results of both obser¬ 
vations, the radius of SN 2011dh, r 2 oiidh, is found to follow 
a power-law, given by r 2 oiidh = 5.85 x 10®® (tage/30days)® ®^ 
cm. At tage = 159 days, r 2 oiidh = 2.71x10®® cm ~ 0.11 mas, 
which is < 0.1% of the synthesised beam of our EVN obser¬ 
vation. 

SN 2011dh has been observed intensively with the 
JVLA across several frequency bands b etween 1.4 to 43 GHz 
jKrauss et alJ[201^ : lHoresh et al.ll2013l l between = 4 to 
tage = 93 days after to- With a detection of SN 2011dh 


at 159 days after to, it may be pos sible to provide con - 
straints to the light curve model of iHoresh et_alj _([20ljL 
Figure 1111 includes the 4.9/5 GHz data from iKrauss et al.l 
(I2OI2I) and iHoresh et all ll2013l L with the 1.65 GHz flux den¬ 
sity values determined via the spectral index b etween the 1.4 
and 1 .8 GHz flux densities at each epoch from iKrauss et al] 
1I2OI2II . The lines are model fits to the light curves, which 
are described below. 

The radio emission from a supernova (SN) arises 
through the interaction of the SN ejecta with th e circum- 
stellar medium (GSM) llWeiler et al.l[l98^ . 1200^ 1. The re¬ 
sulting radio light curve is then a result of the competing 
effects from the declining non-thermal radio emission and 
the more rapid declining thermal (free-free) and non-thermal 
synchrotron self-absorption (SSA) as the ejecta propagates 
through the GS M. This relationship is generalised in a para¬ 
metric model bv lWeiler et al.l (1200211 . which allows for power- 
law variations in key quantities and can be written in the 
simplified form, 


F 


K1 


( V / tagB -to\ 

V5GHz/ V Iday ) 




1 - 


^opt 


mJy, 


(5) 


where a describes the frequency dependency (or spec¬ 
tral index), TLl is a scaling factor and is related to the 
flux density at 5 GHz 1 day after the explosion, and /3 de¬ 
scribes the powe r-law dependency o f the flux density with 
the SN age, t^ge JWeiler et al.l 120021 . The absorption term, 
Topt, describes an internal absorption by material mixed with 
the emitting component, which assumes a planar geome¬ 
try (IWeiler et al.ll200^ : iHoresh et al.ll2013fl . Both the early 
light curve {t < 15 days) and the radio spectra {t < 93 
days) were found to be consistent with abso rption due to 
SSA (iKrauss et al.|[^12l : iHoresh et aLll2013ll . The absorp¬ 
tion mechanism in Equation [5] is therefore taken to be of the 
non-thermal SSA form and is defined by the optical depth, 
Topt, given as. 



where A'5 is the optical depth 1 day after the explosion 
and 5 describes t he power-law depen dency of Topt with the 
SN age, tags — to (IWeiler et al.ll2003l . 

The d ashed lines in Eig 1111 are the fit to Equation [5] 
derived by IHoresh et al.l (l2013l l. using 22.5 GHz JVLA data 
for t < 15 days. While the model from IHoresh et all (l2013fl 
provides a better fit to the observed data for the early rise 
of the light cure, it clearly over-estimates the flux density 
for the later times. Using the 4.9/5 GHz data, a new SSA 
model fit to the light curve was obtained for Equation [SI 
which provides a better overall fit to the data at 1.65 GHz, 
including the turnover observed with the EVN. The final fit 
parameters are a = -1.30, f5 = -1.14, J=-2.88, K1 = 928.5 
mJy, and K5= 5.5 x lO"* and are displayed in Eigure fTT] ('the 
solid lines). 

The difference between the model fits can be explained 
through the estimates of Topt, which are displayed in Fig¬ 
ure 1121 At the early stages of the SN expansion, the GSM 
is optically thick (Topt ^ 1) and gradually becomes opti¬ 
cally thi n (Topt ^ 1) as the SN expands into the surround¬ 
ing ISM (IWeiler et al.ll2003 L coinciding with the emission of 















































































16 H. Rampadarath et al 



Figure 11. Light curve of SN 2011dh at 5 GHz (large filled cir¬ 
cles) and 1.65 GHz (small filled circles and square). The 5 GHz 
flux density was o btained from the JVLA (IKrauss et al.l [20121: 
iHoresh et al.|[2013h . while the 1.65 GHz flux density data are 
from the JVLA (small filled circles) and the EVN observation 
(square) from this study. The lines are least- square fits to Equa- 
tionO from this study (solid lines) and from fHoresh et ^ ll2013l . 
dashed lines). The arrows are upper limits on the flux de nsity (left 
- 5 GHz and right - 1.65 GHz^from lHoresh et al.l <l2013l) obtained 
with the JVLA at these times. 



Figure 12. SSA optic al depth, obta ined via Equation [6] 

using the parameters bv iHoresh et al.l l|2013h . (solid line) and this 
study (dashed line) at 1.65 GHz. 


the peak radio brightness. Furthermore, the time at which 
the medium becomes optically thin is dependent on the fre¬ 
quency of emission . Thus, in deriving the model fit using 
the 22.5 GHz data. [Horesh et qH l|2013l) underestimated the 
optical depth (Figure [T^ blue dashed lines) and, therefore, 
the attenuation of the flux density at the lower frequencies. 

X-Rays. SN 2011dh was clearly detected for all obser¬ 
vations following the explosion date. Similar to SN 19941, 
the flux of SN 2011dh was found to be decreasing with time. 
T he Chandra d a ta lis ted in Table [3] were used extensively 
bv iMaeda et al.l ll2014l l to investigate the long-lasting X-ray 
emission from SN 2011dh and the mass loss history of the 
yellow supergiant progenitor. We find that our flux measure¬ 
ments are consistent with their reported values. 



13'’30'"l 1.014^ 11,013^ 11,010’' 11,008' 11,006' 

Right Ascension (J3000) 

Figure 13. A higher resolution EVN image of 30114-1041 made 
with Briggs robustness = 0, cellsize of 0.1 mas. Beam size = 
6.23 mas X 3.19 mas, position angle = 13.7° and image RMS = 
21.9 /rjybeam~^. 

5.4 Background Radio Sources 

54.1 JI33OII+47IO4I 

The radio source J133011-I-471041 was detected with the 
VLA at both 20 c m (500T22 /rJy beam"^) and 6 cm (482±50 
pjybeam“^) bv iMaddox et aU l|2007l l. giving a flat radio 
spectrum {a = -0.03). Comparing the radio position with 
archived Chandra X-ray data, they found a 10®^ ergs s”*^ 
counterpart, making it the brightest X-ray/radio overlap 
in their survey. However, no Ha emission was found at 
the source po s ition from the HST ACS archive, leading 
iMaddox et al.l ll2007l l to suggest that J1330114-471041 may 
be a microquasar in a radio-loud state during the time of the 
VLA observations. If indeed this source is a microquasar, its 
luminosity at 5 GHz would be a factor of two more luminous 
than the most luminous discovered to date (G RO J1655-40: 
iHiellming fc Rup^ll995l . iTingav et al.lll99^ . 

The VLBI image (Figure [^) reveals a compact, rela¬ 
tively unresolved source with a total EVN flux density (see 
Table[T]) matching the VLA flux density within measurement 
errors. Since there is little or no extended emission resolved 
by our VLBI observations, the radio emission seen with the 
VLA and EVN is emanating from the same compact region, 
a few mas in size. Re-imaging the VLBI data with a higher 
resolution by setting the Briggs robust parameter to 0 and 
adopting a cellsize of 0.1 mas, reveals a double morphology, 
resembling a core-jet structure typically seen in radio galax¬ 
ies or AGNs. The resulting image is shown in Eigure 1131 
with the peak surface brightness of the components of 270 
pjy beam“^ and 233 ^Jybeam“'^, separated by 4.3 mas. To 
the east of the double, at a distance of 38.6 mas from the 
brightest component of 3011+1041, is a weak unresolved 
source with peak surface brightness of 123 /ijybeam“^ that 
may be a jet hotspot related to J133011+471041. 

X-ray Spectrum. The stacked Chandra spectrum (855 
ks) is well fitted with a simple power-law (Eigure I14|l with 
photon indeL^ T = 1.68lQjg, and line-of-sight Galactic ab¬ 
sorption Nh — 2x 10^° cm“^. Thermal models such as disk 
black-body or thermal plasma are ruled out; adding a soft 
thermal component to the power-law model does not im¬ 
prove the fit. The observed flux of the stacked spectrum is 
/ 0 . 3-10 = (2.4 ± 0.1) X 10“^^^ erg s“^ cm“^, from which we 

r is related to the spectral index a commonly used in radio 
astronomy by F = 1 + a, where S oc v~°‘. 
























































A High Resolution Wide-Field Radio Survey of M51 17 



0.5 1 2 5 

Energy (keV) 


Figure 14. Stacked Chandra/ ACIS-S X-ray spectrum of 
.1133011+471041 fitted with a power-law model (F = 1.7 ± 0.1, 
and line-of-sight absorption), and corresponding residuals 
(xl = 100.23/105 = 0.95). 


infer an emitted luminosity Lo, 3 _io = (2.1 ± 0.1) x 10^® erg 
if the source belongs to M51a. In fact, we suggest that 
this source is more likely to be a background AGN. 

We can estimate the likelihood that J133011+471041 is 
a background source using the fundam ental-plane relation 
dMerloni et al.ll2b03l : IPlotkin et al.ll2012l ). If it were in M51a, 
logLR(5GHz) « 34.2 a nd logZ/ 2 -ro ~ 38.1. Using the best¬ 
fitting coefficients from l^otkin et ^ ll2012tl . the resulting 
black hole mass would be ~ 10® Mq , inconsistent with an X- 
ray binary or even an intermediate-mass black hole in M51a. 
The only possibility for J133011+471041 to be located in 
M51a is if its radio flux is strongly Doppler boosted, which 
we cannot yet rule out. 


Optical/X-Ray. The ratio between X-ray and opti¬ 
cal luminosity is consider ed to be a relia ble indicator of X- 
ray source classif ication IIGo mastril 120041). The relationship 
is defined as (se e Maccacaro et al.lll988l : IBi 

iMcHardv et al.ll2003l ): 


Tsarger et al.ll 200 ' 2 : 


logio(-Px/-F’opt) = logioTx -I- 5.5 -I- Rmag/2.5, (7) 

where Fapt is the optical flux (in cgs) and J?mag is the 
Johnson-Gousins R-band flux density. The distribution of 
login(Tx/ToDt ) for spectroscopicall y identifi ed X-ray AGN 
from ROSAT (|Hasinge£_et_alj 19981) . A SCA llAkivama et akl 
l2003l). G/tandra jGiacconi et al.l 2001J), and XMM-Newton 


iMainieri et al. I I2OO2I: iFiore et al 


20031) surveys fall within 


-1< logjQ(fx/f’opt) < 1, while for stellar mass sources, such 
as x-ray-binaries and microquasars, logpQ{Fx/Fapt) > 2. 

Since the galaxy was not observed in the R-band with 
the HST/AGS, to estimate the i?mag of J133011+471041, 
the fluxes of the HST B, V and I bands were first estimated 
as described in Sectionjd] The R-band flux (10.26 /iJy) was 
then estimated via a non-linear least squares fit to the flux 
densities of the B, V and I bands, from which a value of 
21.1 for Rmag was determined. Substituting the values into 
Equation[7]results in log,^Q(Tx/T’opt) = -0.03, indicating that 
3011+1041 is likely a background AGN and is consistent 
with what was found using the fundamental-plane argument 
and the X-ray spectrum. 

Radio Loudness and AGN Classification. Since 
we have concluded that J133011+471041 is a background 


source, based upon its optical/radio/X-ray flux ratios, we 
can now classify its AGN type from its radio loudness para m- 
eter Rl, defined as llTerashima fc Wilsorj|200^ : lHoll2008l ): 


Rl ~ 


vFR5 GHz) 
/x(2-10keV)' 


( 8 ) 


Substituting f/(5 GHz) = 482 pjy and / 2 _io = 1.5 x 
lO”^'^ erg s“^ cm“^ gives log i?i, « —2.8. There¬ 

fore, J133011+471041 falls within the parameter space 
of radio-loud low-lumino sity AGNs (se e Figure 4 of 
iTerashima fc WilsorJl2003l and Figure 10 of lHd[200^ ). 


5 . 4.2 J133005+471035 

The radio source J133005+471035 is resolved into four com¬ 
ponents with the FVN at 18 cm (Figure[6]3), where the com¬ 
ponents are labelled G1 to G4 in order of decreasing peak 
surface brightness. 

Table [5] lists the different properties of the four compo¬ 
nents, including the position of the peak surface brightness, 
peak and total flux densities, size of the component’s ma¬ 
jor axis and brightness temperature. The VLBI morphology 
of J133005+471035 shows evidence for a two-sided jet ema¬ 
nating from an active nucleus. However, without higher fre¬ 
quency VLBI observations it is difficult to determine which 
is the core that hosts the active nucleus. The similarity of 
both G1 and C2, in terms of size, flux density and bright¬ 
ness temperature suggests that they are, most probably, ra¬ 
dio lobes/hotspots resulting from jets emanating from an 
unresolved, fainter active nucleus. 

In lower resolution radio maps of M51a with the VLA 
and MERLIN, J133005+471035 is the brightest radio source 
in the M51 system. Figure [15] shows the VLA-A 20 cm and 
MERLIN 18 cm maps in relation to the EVN 18 cm map 
(represented by the green crosses). In the FIRST survey the 
source is unresolved at 20 cm with a total flux density of 
12.64 ± 0.14 mjy, while at higher resolution VLA-A, a total 
flux density of 13.38 ± 1.50 mJy is recovered from the un¬ 
resolved 20 cm VLA-A image in Fig ure 1151 This is simila r 
to the peak flux density obtained bv iMaddox et all ll2007t) : 
9.6 mjybeam“^ at 20 cm and 4.3 mjybeam“^ at 6 cm, 
giving an intermediate spectral index of —0.66. The MER¬ 
LIN 18 cm image displayed in Figure [15] shows an elongated 
structure, closely resembling the one-sided c ore-jet morphol¬ 
ogy of AGNs and quasars on kpc-scales (]Bridle Sz Perlevl 
1 19841) . The total flux density as measured (via blobcat) 
from the MERLIN image is 9.95 ± 1.00 mJy, recovering 
~ 80% of the flux density of the VLA-A at 18 cm (after 
correcting the 20 cm flux density with a spectral index of 
-0.66). Subtracting the total VLBI flux density listed in Ta¬ 
ble [^ from the MERLIN 18 cm flux density gives 8.12 ± 
1.01 mJy, which may result from thermal and non-thermal 
emission from star formation and non-thermal emission from 
kpc-scale jets. 

A faint X-ray sonrce was found coincident with the posi¬ 
tion of J133005+471035 in the stacked Chandra image (855 
ks). We measure a 0.3-10 keV count rate of 6.0+1.2 x 10~® ct 
s“^, mostly distributed around ~1 keV. Using Cash statis¬ 
tics, we obtain that the spectrum is best fitted by optically- 
thin thermal plasma emission, at kT = O.L/q ® keV, ab¬ 
sorbed by a column density Na « 5 x 10^*^ cm“^ (probably 
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MERLIN 18 cm 



Right Ascension {J2000) 

Figure 15. Radio images of the radio source J133005+471035. VLA-A: Contours begin at ±10% the peak flux density (11.8 mjy beam“^), 
and increases in increments of 10% of the peak flux density. MERLIN: Image was made with the full MERLIN array (i.e. including 
the long baselines to Cambridge), which resulted in a beam size of 0.196” X 0.164” . The contours begin at ±3crinjagg and increase with 
multiples of \/2i where aiYna,ge = 60.3 /ijybeam“^. EVN: The image is a zoomed in version of the one shown in Figure|6l), where Cl - 
C4 are the different components, labelled in order of decreasing peak surface brightness. 


Table 5. VLBI detected components of J133005±471035. Columns: (1) VLBI component defined in Figure [TSl (2) & (3) VLBI source 
position; (4) Angular separation of the Fp of each component from the Fp of Cl; (5) Peak flux density; (6) Integrated flux density; (7) 
Size of the major axis from imfit; (8) Brightness temperature at u MHz derived from the Rayleigh-Jeans formalism; {c^ Fj) / {2 tt , 

where c is the speed of light and k the Boltzmann constant. 


Position of Fp 


VLBI 

Source 

(1) 

RA 

(2) 

Dec 

(47° 10') 

(3) 

©Cl 

(mas) 

(4) 

Fp 

(^Jy beam~^) 

(5) 

Fi 

(hJy) 

(6) 

Lmaj 

(mas) 

(7) 

logic Tb 
(K) 

(8) 

Cl 

5.1381(2) 

35.7896(4) 

0 

211.7T25.6 

457.2T65.2 

12.8 

6.26 

C2 

5.1373(3) 

35.787308 

9.72 

199.8T24.6 

414.0T59.1 

10.98 

6.34 

C3 

5.1355557 

35.782491 

27.5 

99.1T17.3 

313.0T65.3 

14.65 

5.98 

C4 

5.1404641 

35.794635 

24.2 

91.8T16.6 

650.2T128.6 

32.03 

5.61 


due to the location of the source in a dust-rich spiral arm). 
A power-law model is ruled out because it would require an 
unphysically steep photon index F > 5. The observed flux 
corresponding to the best-fitting thermal-plasma model is 
fo.3-10 = 2.llg g X 10“^® erg s~^ cm“^, and the unabsorbed 
luminosity is Lo.s-io = O.SlJ g x 10®^ erg s“^ if located in 
M51a. The thermal nature of the source makes it consistent 
either with an SNR, if the source is in M51a, or with knots 
of hot gas shock-ionized by a jet. 


5.4.3 JI33OI6+47IO24 

VLA-A 20 cm observations of the radio source 
J133016T471024 shows a morphology (see Figure 
similar to near by powerful FR I I rad io galaxies such 
as Cygnus A (ICarilli fc Bartli^ Il996l) and Pictor A 
dPerlev et al.l Il997l) , displaying the morphological compo¬ 
nents common to this type of galaxies. These sources are 
composed of a compact, high surface brightness temper¬ 
ature, typically flat spectrum source at the centre, from 
which elo ngated radio jets end in high surface brightness 
hotspots llCarilli &: Bart hell 1 19961 ). The source was detected 
as a partially resolved double radio source in the FIRST 
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VIA FIRST VLA-A 20 cm 



Figure 16. Radio images of the source J133016+471024. VLA: the contours of both the FIRST and VLA-A images begin at ±10% the 
peak flux density (6.02 mjybeam~^&; 0.982 mjybeam”^, respectively), and increases in increments of 10% of the peak flux density. 
MERLIN: the first contours are at ±3 it, and increases in multiples of \/2, where a = 50/rJy beam“'^; EVN: a zoomed in version of the 
one shown in Figure |6] 


survey with total flux densities of 12.56±0.63 mjy and 
11.67±0.58 mJy for the southern and northern lobes, 
respectively. The higher resolution VLA-A 20 cm image 
resolves the double sources into two radio lobes and a 
central radio source, with evidence for elongated radio jets 
ending in hotspots at the extremities of the radio lobes. The 
VLA-A image recovered 82% and 71% of the FIRST flux 
density of the southern and northern lobes, respectively. 
The central source was found to have a total flux density of 
2.15±0.04 mJy at 1.4 a-nd 0.98±0.01 mJy at 5 GHz, 

giving a spectral index of -0.62. 

The central radio source is unresolved with the higher 
resolution MERLIN 18 cm observations, with total flux den- 


The 1.4 GHz dataset was re-imaged with the same beamsize 
as the 5 GHz dataset. 


sity of 0.97T0.14 mJy. However, there is evidence for low sur¬ 
face brightness emission (~150 /rjy beam“^) along the direc¬ 
tion of the southern lobe/jet (see Figure fTBll . and at the lo¬ 
cation of the peak VLA-20 cm surface brightness component 
of the northern lobe (not shown in Figure [161) . The VLBI 
observation only detects the central radio source, showing 
a slightly elongated structure extending from the compact 
source in the direction of the southern jet, suggesting that 
the southern jet/lobe may be relativistically beamed to¬ 
wards our line of sight. From the total flux density and 
major axis size of the central source we obtain a bright- 
ness temperature of logi p Tb = 7.1, typical o f AGNs (e.g. 
iBlundell fc Beaslevlll998 : lUlvestad et alJl2005fl . 

In the standard FR I/FR H classification model 
llFanaroff fc Rile^ll974l ) both types of radio sources can be 
divided by radio luminosity, where the dividing luminosity 
is 10^® W Hz“^ at 1.4 GHz (below lies the FR I and above 
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Figure 17. Stacked Chandra/ACIS-S X-ray spectrum of 
3016-1-1024 fitted with a power-law model (P = 2.0di0.2, and line- 
of-sight a bsorption). The unbinned data were fitted with Cash 
statistics llCas ainii), and they were then rebinned to a mini¬ 
mum signal-to-noise ratio of 2.5 for display purposes only. 

the FR II galaxies). Based upon the morphology shown 
in Figure [TSl it is possible to classify J133016+471024 as 
an FR II galaxy. Assuming a minimum luminosity of 10^® 
W and with a total flux density at 1.4 GHz of 2.15 

mjy, J133016+471024 would be at a minimum distance of ~ 
6246 Mpc {z ^ 0.95). Note, placing Cygnus A at 6246 Mp c 
would result in a total flux density ~5 mJy at 1.4 GHqfj. 

From the Sloan Digital Sky Survey (SDSS) Data Re¬ 
lease 9 llAhn et al.ll2012ll . an optical source was found to be 
coincident with the position of the central radio source of 
J133016+471024. The R-band optical data were obtained 
through the SDSS Science Archive Server in the “corrected 
frame" format, which have been calibrated in nanomaggiej^ 
per pixel, and have had a sky-subtraction applied. We per¬ 
formed aperture photometry with standard Heasoft pack¬ 
ages and found a total flux density of 13.11 nanomaggies or 
47.6 /rjy is obtained, which is equivalent to an Rmag of 19.5. 

Inspecting the stacked Chandra dataset, we found an 
X-ray source at the position of the central radio source of 
J133016+471024. We extracted a combined spectrum as de¬ 
scribed in Sect l4.3l The spectrum is best fitted (Figure [TT]) 
with a simple power-law of photon index F = 1.98lQ jg 
and line-of-sight Galactic absorption Ah = 2 x 10^'^ cm“^; 
adding free intrinsic absorption does not improve the fit. 
The observed flux is / 0 . 3-10 = 4.2^q'4 x 10“^® erg s“^ 
cm“^. Using Equation O we determine a radio loudness pa¬ 
rameter Rl = —1.97, which places this background source 
firmly in the regime of r adio-loud AGN ( see Figur e 4 of 
iTerashima fc Wilsonll200lj and Figure 10 of lHollioO^ . 

5 . 4.4 J132932+471123 

While J132932+471123 is not the weakest VLBI source in 
our sample, it has the lowest S/N ratio (see Table [T]) which 
results from a high RMS noise due the source’s position 

Assuming the Flux densi ty of Cygnus A is F i. 4 GHz = 1598 Jy 
at a distance of 247.2 Mpc llBirzan et ahll^OlTl . 

A star of brightness 1 nanomaggie has a magnitude 
of 22.5 in any band, or a flux density of 3.631 pJy. 
(WWW. sdss3. org/dr9/help/glossary, php) 


at the edge of the primary beam of Ef. Its positional off¬ 
set is 3 ' .49 in RA and 19” .6 in Dec from the M51a cen¬ 
tre, placing it at the edge of th e M51 a disk and beyond 
the survey area of iMaddox et al.] (l2007l ). At the VLBI po¬ 
sition of J132932+471123, a radio source is detected in the 
VLA-A 20 cm image with a peak flux density of 237+49 
/rJybeam“^. However, no corresponding source was found 
in the VLA-B 6 cm image above 60 /ijybeam~^. This 
suggests that J132932+471123 is a steep spectrum source 
(a < —1.10, S oc v°‘). 

In all the Chandra observations, the source is located a 
few arcmin from the aimpoint, near the edge of the S3 chip, 
where the point spread function is larger (several arcsec) and 
it becomes more difficult to detect faint sources. Nonethe¬ 
less, we find a 3 (t detection in the combined X-ray image 
(835 ks), with a 0.3-10 keV count rate of (3.2 ± 1.4) x 10~® 
ct s“^. A Gash-statistic power-law fit gives a photon index 
F = 2.0lj g with fixed line-of-sight absorption; there are not 
enough counts for any additional fit components. The ob¬ 
served flux / 0 . 3-10 = (5 ± 2) X 10“^® erg s“^ cm~^, which 
corresponds to an X-ray luminosity « 5 x 10®® erg s“® if 
the source is located in M51a. No optical counterparts were 
found in any of the archival images used for this study; this 
is not inconsistent with an AGN identification, because we 
expect R « 24.5 mag for an X-ray-to-optical flux ratio of 1. 
In summary, considering its location at the outskirts of the 
galactic disk, we tentatively classify this object as a back¬ 
ground AGN. 

5.5 Estimating the Supernova and Star Formation 
rates 

5.5 .1 Global supernova and star formation rates 

The supernovae rate, vsn of the entire galaxy can be esti¬ 
mated using the age and total number of SNRs in M51. Both 
estimates are given in Sect. 15.H where the ages of the SNRs 
were estimated to be 280 - 950 yrs. As discussed in Sect. 15.H 
our VLBI observatons detected only 3/32 radio sources in 
M51. Since the 3 detected radio sources are background 
AGNs and the AGN nucleus of M51, this suggests that 
the 29 undetected radio sources were SNRs that have ex¬ 
panded beyond the maximum size detectable with our obser¬ 
vations. This suggests a supernovae rate range of 0.03 - 0.10 
yr“®. The global massive star-formation rate (M ^ 5Mq) 
can be estima ted from the vsn using Equation 20 from 
ICondonI lll992lL Thus we obtain SFR{M ^ 5 Mq) ~ 0.73 — 
2.43 Moyr"® {SFR{M ^ 0.1 Mq) ~ 3.65 - 12.15 Mgyr"®). 
The lower limit of our estimation is slightly higher than 
the SFR the total star formation rate of the entire galaxy 
(3.4 Mqyi~^) derived via infrared and ultraviolet observa¬ 
tions llCalzetti et al.ll2005l L This suggests that some of the 
29 undetected radio sources are not SNRs. A more robust 
estimate on osn and the SFR from the radio observations 
can be achieved followi ng the multi-epoch, multi-f requency 
Monte-Garlo method of iRampadarath et al.l ll2014ll . requir¬ 
ing further observations. 

5. 5.2 Star formation rate of the inner nucleus 

The star formation rate (SFR) of a galaxy is directly propor¬ 
tional to its radio luminosity, Li, at wavelength u, assuming 
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that all the radio emission is associated with star forma¬ 
tion llCondonI Il992l : iHaarsma et al. IH). An estimate of 
the flux density associated with star formation within the 
inner nucleus of M51 can be obtained by subtracting the 
compact VLBI components from the MERLIN data. A to¬ 
tal flux density of 2.68 ± 0.14 mjy was recovered from the 
MERLIN 18 cm observation, and 0.25 ± 0.02 mJy from 
the EVN observation. Subtraction yields a flux density for 
the extended radio emission at 18 cm of 2.43 ± 0.14 mJy. 
However, there may be contributions to the extended emis¬ 
sion from the jet components that were resolved-out by the 
EVN observations. Since this cannot be distinguished by 
the MERLIN observations, the flux density of the extended 
emission places an upper limit of the radio emission associ¬ 
ated with star formation and hence the SFR _ 

Using the stand ard r elations given i n CondonI lll992tl . 
IHaarsma et al.l (|2000| L and|CMndpn_et^D |l20 p^, along with 
the implied standard IMFs I Saloeteil Il955l : iMiller fc Scald 
[l9Zi), we derive an upper limit for the star formation 
rate of SFR(M ^ 5 Mq) < 0.044 Mgyr-^ (or SFR{M > 
0.1 Mq) < O.24M0yr“^) for the inner nuclear region of 
M51, which is ~10% the total star formation rate of the 
entire galaxy. 


6 CONCLUSIONS 

This paper presents the deepest, widest VLBI survey of a 
nearby grand design spiral galaxy. The target for our study 
is the Whirlpool galaxy (M51a, NGC 5194), that is un¬ 
dergoing a merger with its smaller companion, NGC 5195. 
The full disk of M51 was surveyed using the multi-phase 
centre technique (192 phase centres) for compact radio 
sources. By combining this single pointing of the EVN, 
with multi-wavelength radio observations. X-ray and opti¬ 
cal data we investigated: the properties of the nuclear re¬ 
gion of M51a (including a pc-scaled jet); the light curves of 
the Type Hb supernovae, SN 201 Idh (radio) and the Type 
Ib/c supernova, SN 19941 (radio and x-ray); the nature of 
four background AGNs (J133005+471035, J133011+471041, 
J132932+471123, J133016+471024); and the supernova and 
star formation rates of M51a. 

It is important to note that of the 192 phase centres, 
only six provided positive detections, indicating that most of 
the computing time was spent on imaging empty fields. Fu¬ 
ture VLBI surveys of nearby galaxies would therefore benefit 
from limiting the phase centres to only sources detected with 
lower resolution interferometers (e.g. the JVLA and eMER- 
LIN). Nonetheless, this study demonstrates the extensive 
scientific capabilities of a single wide-held VLBI observa¬ 
tion. 
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APPENDIX A: LR METHOD 

The ratio used in Sect. l3.2T^ to estimate the likelihood that a 
pixel in a VLBI image is associated with a catalogued source 
(the alternative hypothesis, Hi), or a noise spike (the null 
hypothesis, Hq) is given as, 


Table Al. Parameters used in the likelihood ratio method to 
search for sources within the EVN images. Columns: (1) Cata¬ 
logue of radio sources in M51; (2) Number of sources detected in 
M51 by both catalogues; (3) Major axis of the survey beam; (4) 
Number of pixels per dimension required to image the field of size 
Bmaj • 


LR = 


p(gi) 

p{H^) ’ 


(Al) 


where p(-ffo) is the likelihood the pixel is a noise spike 
and is taken as the Cumulative Distribution Function of the 
positive half of the Gaussian distribution, while p{Hi) is the 
likelihood the pixel is associated with a catalogued source 
and is dependent on the position of the pixel in relation 
to a catalogued sou rce position/morphology as defined by 
[Morgan et alj ll2013l l. Before applying the LR method, the 
AlPS task UVFIX was used to produce visibilities centred on 
the catalogued positions allowing the data to be averaged 
in frequency and time, with little loss of amplitude response 
due to image smearing. Dirty images were generated for each 
of the catalogued source positions with the full arra-yPI. using 
parameters: cellsize = 0.5 mas; image size listed in Table I^Tl 
and natural weighting. In addition, the outer 512 pixels re¬ 
moved prior to applying the LR method. The required VLBI 
parameters to search for sources from both catalogues are 
listed in Table 

LR images were generated for each catalogued source 
and the position of the maximum LR value was visually 
compared to the corresponding VLBI peak pixel S/N. Pix¬ 
els with logic (LR) > 10 and S/N > 6.7a were taken as 
confirmed detections. Tables IA2I and IA3I list the LR values 
for the confirmed detections along with the corresponding 
VLBI image parameters for confirmed detections, and the 
non-detection with t he highest LR value for both surveys. 
As with the results of iMorgan et al.l ll2013tl . the LR test was 
found to be more dependent on the VLBI pixel S/N (i.e. 
p{Ho)), than the positional constraints given by the lower 
resolution catalogues. 

An example of this is displayed in Figure lAll where 
a 22a VLBI pixel was discovered from two FIRST 
sources, but has logic (LR) values >100 when compared 
with both FIRST sources. However, pixels that were found 
to have better positional constraints (i.e. higher values of 
nevertheless, resulted in lower LR values. As such 
no significant advantage is given by the LR method over 
simply detecting the brightest pixels and thus, it is difficult 
to determine whether pixels with S/N< 6.7a are indeed real 
detections, and are henceforth considered as resulting from 
noise. _ _ 

Of the three [Maddox et all (l2007l l sources detected by 
the EVN, two were found with very good agreement be¬ 
tween the VLBI and lower resolution source positions (see 
Figure lATb L It ma y be possible to use e stimates of the po¬ 
sitional error of the lMaddox et all ll2007t l survey to improve 
the LR method and increase the chance of detecting sources 
below 6.7a. However, this is beyond the scope of this study. 



Source 

Bmaj 

Image size 

Catalogue 

Count 

(") 

(pixels) 

(1) 

(2) 

(3) 

(4) 

FIRST . 

15 

5 

7936 

Maddox (20 cm) 

96 

1.51 

1984 


This test was also conducted for the reduced array of maxi¬ 
mum baseline of 7.5 MA. However there was no difference in the 
final result compared to the full array. 
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Table A2. Results of the likelihood ratio analysis of the VLA First catalogue. Columns: (1) Flux density of the brightest pixel from 
images made with the full EVN array; (2) The Icr of the individual images obtained through fitting a Gaussian distribution to the 
histogram of pixel fluxes; (3) Signal to noise ratio; (4) The probability the brightest pixel in the image is a noise spike, i.e. a source of 
S/N ratio would occur by chance once in every l/10~^°9io{p{^o)) pixels; (5) Likelihood the pixel is related to a catalogued source based 
upon position; (6) The likelihood ratio parameter of the brightest pixel in the image. See text for details. 


VLBIpeak 
(/4jy beam“^) 
(1) 

cr 

(/.tjy beam“^) 
(2) 

S/N 

(VLBIpeak/^) 

(3) 

-logio(p(Ho)) 

(4) 

-logic (p(Hi)) 

(5) 

logic (LR) 

(6) 

326.42“ 

11.75 

27.77 

169.4 

0.49 

168.9 

233.00*’ 

10.03 

23.24 

118.9 

0.89 

118.0 

231.16*’ 

10.09 

22.90 

115.7 

1.09 

114.6 

84.13“ 

8.96 

9.39 

20.5 

1.09 

19.4 

53.35t 

9.00 

5.93 

8.8 

0.60 

8.2 


Source identification (see Table ITl: “ J123005+471035; ^ A single VLBI source (J133016+471024) located between two FIRST 
sources. See Figure fAll ^ J132952+471142. ^ The non-detection with the highest LR. 

(6) The likelihood ratio parameter of the brightest pixel in the image. See text for details. 


Table A3. Results of the likelihood ratio analysis of the Maddox catalogue. 


VLBIpeak 
(/ijy beam“*^) 
(1) 

(J 

(/ijy bearn”^) 
(2) 

S/N 

(VLBIpeak/<7) 

(3) 

-logio(p(Ho)) 

(4) 

-logio(p(Hi)) 

(5) 

logic (LR) 
(6) 

327.26“ 

17.33 

18.31 

74.5 

0.49 

74.0 

317.42*’ 

30.23 

10.82 

26.9 

0.40 

26.5 

86.07“ 

10.12 

8.51 

17.0 

0.30 

16.7 

47.48t 

8.52 

5.57 

7.9 

0.60 

7.3 


Notes. Symbols and headings have same meaning as Table rA2| 

Source identification (see Table[T}: “ J133005+471035; ^ J133011+471041; J132952+471142. 



Figure Al. (a) The location of the peak VLBI pixel with respect to the catalogued position of the nearest corresponding FIRST source. The circles are: 
red (J133005-|-471035); and blue (J132952+471142). (b) The location of the single VLBI source (J133016-1-471024; white cross) located between two FIRST 
sources (greyscale). The ellipse gives the FIRST beam size (see Table fAH . (c) The location of the peak VLBI pixel with respect to the catalogued position of 
the nearest corresponding 20 cm iMaddox et al ] i2007l 1 source. The circles are: red (J133005+471035)i green (J133011+471041); blue (J132952+471142). The 
ellipse in both (a) and (c) are the beam sizes of the respective surveys, centred on the catalogued position. The size of the circles indicate the astrometric error 
on the catalogued positions via the relationship, • crl2Sp, where Bmaj/min major or minor beamsize, cr is the RMS noise in the catalogue 

images and Sp is the peak flux density of the catalogued source ^Morgan et alJ|2013D. Dashed circles have been placed around the sources to aid the reader. 





















